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ABSTRACT 

We have conducted a survey of a sample of infrared-dark clouds (IRDCs) with the Spitzer Space Telescope 
in order to explore their mass distribution. We present a method for tracing mass using dust absorption against 
the bright Galactic background at 8 /vm. The IRDCs in this sample are comprised of tens of clumps, ranging 
in sizes from 0.02 to 0.3 pc in diameter and masses from 0.5 to a few IQ^Mq, the broadest dynamic range in 
any clump mass spectrum study to date. Structure with this range in scales confirms that IRDCs are the the 
precursors to stellar clusters in an early phase of fragmentation. Young stars are distributed in the vicinity of 
the IRDCs, but the clumps are typically not associated with stars and appear pre-stellar in nature. We find an 
IRDC clump mass spectrum with a slope of q'=L76±0.05 for masses from 3QMq to 3000Mq. This slope is 
consistent with numerous studies, culled from a variety of observational techniques, of massive star formation 
regions and is close to the mass function of Galactic stellar clusters and star clusters in other galaxies. We assert 
that the shape of the mass function is an intrinsic and universal feature of massive star formation regions, that 
are the birth sites of stellar clusters. As these clouds evolve and their constituent clumps fragment, the mass 
spectrum will steepen and eventually assume the form of the core mass function that is observed locally. 

Subject headings: ISM: clouds, dust, extinction — stars: formation 



1. INTRODUCTION 

Our understanding of star formation has grown primar- 
ily from the study of local regions forming low-mass 
stars in relative isolation. Nearby regions suc h as Taurus 
(e.g. [Goldsmith etani2008h. p Ophiuchus (e.g. Young et al 
200m~ Johnstone et alj~ 20001). Perseus (e.g . Enoch et al . 
2008"; 'J0rgensen et al."2007'; 'Kirk et al."2006'). Serpens (e.g. 
lesti & Sargent 1998; Harvev et al. 2006; Enoch et al. 200%^ 
the P ipe nebula (e.g. Lombardi et al. 2006; Muench et al] 
12001) . and Orion (e.g. Li et aL.2007: Johnstone et aL 200ll) 
have been studied to great lengths using a variety of tech- 
niques including dust emission, extinction mapping, and 
molecular line emission. This past decade of research has 
shown that star formation regions are assembled hierarchi- 
cally. Within molecular clouds (tens of parsecs in size, 
co ntaining 10"^- 10^ Mp), w e adopt the nomenclature used 
by iBergin & Tafallal (12007^ distinguishing "clouds" (10^ - 
lO"* Mo, 10° - 10^ pc), "clumps" (10-10^^ Mo, lO^'-lO" pc), 
and "cores" (10"'-10' Mo, lO^^-lO"' pc). 

In studies of nearby regions, it is possible to re- 
solve pre-stellar core s with single dish observations (e.g. 
iJohnstone et alj|2000 !). This permits the examination of the 
properties of the fragmentation of the natal molecular clouds 
into smaller components. It is then straightforward to con- 
struct a mass function of cores and a mass function for indi- 
vidual cores can then be constructed. The core mass distribu- 
tions typically derived from dust emission studies are found 
to be strikingly similar to the mass spectra of stars, implying 
that the masses of stars are a direct result of the way in which 
the natal molecular cloud fragments. In contr ast, when CO 
line e mission is used as a mass probe for cores (iKramer et alJ 
Il998h . a more top-heavy distribution results. 

While these studies have brought us a deep understanding 
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of isolated low-mass star formation, this is not the complete 
picture for star formation in the Galaxy. The necessary ingre- 
dient for star formation is dense molecular gas. However, the 
H2 distribution in the Milky Way is not uniform. The primary 
reservoir is in the Molecular Ring, which resides at 4 kpc from 
the Galactic Center a nd contains ~7 0% of the molecular gas 
inside the solar circle dJackson et al.u2006) . Thus, the Molec- 
ular Ring is the heart of Galactic star formation. Indeed, as 
iRobinson et alJ (1 19841) show, the peak of Galactic far-infrared 
emission originates from this region. 

In local clouds, most of the recent progress in our un- 
derstanding of the early stages of star formation have come 
from studies of pre-stellar objects. Populations of starless 
cores have been identified in numerous local regions, and 
they are uni versally cold and quiescent, exhibiting thermal 
line widths (Bergin & Tafalla 2007). Recently, a population 
of cold, dense molecular clouds within the Molecular Ring 
were detected against the brig ht Galactic mid-infrared back - 
ground (from 7 to 25 fim; Egan et alJl998l:ICarey et alJl998h . 
These clouds are opaque to mid-infrared radiation and show 
little or no typical signs of star formation, such as association 
with IRAS point sources. Initial studies demonstrated that 
these objects, termed infrared-dark clouds (IRDCs) are dense 
(n(H2) > 10^ cm-3), cold (T < 20K) concentrations of 10^ 
- 10^ Mo of molecular gas. Based upon the available mass 
for star formation, infrared-dark clouds are likely the sites of 
massive star formation. 

Since their discovery, further studies of infrared-dark 
clouds have established their place as the precursors to clus- 
ters. A number of studies have detected the presence of 
deeply embedded massive protostars using sub-millimeter 
probes (Be uther & Steinacker 2007; Rathb orne et alj |2007|; 
iPillai et al.li2006bh . which confirms that IRDCs are the birth- 
sites of massive stars. Detailed molecular surveys show that 
molecules su ch as NH3 and N2H^ trace the den se gas ex- 
tremely well (iRagan et al.l 2006; Pillai et al."2006a). as seen in 
local dense prestellar cores (Bergin et al. 2002). Furthermore, 
the molecular emission corresponding to the absorbing struc- 
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ture of infrared-dark clouds universally exhibit non-thermal 
linewidths on par with massive star formation r egions. Other 
studies hav e uncovered t he presence of masers ("Beuth er et al. 
|2p02; Wan g et al.ll2007l) and outflows ( Beuther & Sri dharari 
|2007), known indicators of ongoing embedded star formation. 
Already, the evidence shows that these are the sites where 
massive stars and star clusters will form or are already form- 
ing. In order to understand massive star formation, and thus 
Galactic star formation, it is crucial to understand the struc- 
ture and evolution of IRDCs. 

Studies of infrared-dark clouds to date have left the funda- 
mental properties of cloud frag mentation go relatively unex- 
plored. iRathborne et al.1 (12006) showed that IRDCs exhibit 
structure with median size of ~0.5 pc, but observations of 
IRDCs with the Spitzer Space Telescope, which we describe 
in ^ reveal that there exists structure well below this level. 
We characterize the environment in ^and highly structured 
nature of infrared-dark clouds in ^ by utilizing the high- 
resolution imaging capabilities of the Spitzer. In ^ we an- 
alyze the IRDC absorbing structure, derive the clump mass 
function, and put the results in the context of previous stud- 
ies. We find the mass funct ion to be shallo wer than Salpeter 
initial mass function (IMF) (ISalpeterlll955h and more closely 
aligned with that observed using CO in massive star form- 
ing regions. Given the strong evidence for fragmentation and 
star formation characteristics of these objects, we suggest they 
are in the initial stages of fragmentation. The conclusions as 
well as the broad impact of these results are discussed in ^ 
The results of this study provide an important foundation for 
further studies of IRDCs with the instruments of the future, 
allowing us to probe the dominant mode of star formation in 
the Galaxy, which may be fundamentally different from the 
processes that govern local star formation. 

2. OBSERVATIONS & DATA REDUCTION 
2.1. Targets 

Searc hing in the vicinity of ultr a-compact HII (UCHII) 
regions ("Wood & Churchweiri989^ for infrai-ed-dark cloud 
candidates, Ragan etal. (2006) performed a survey of 114 
candidates in N2H+(l-0), CS(2-1), and C"*O(l-0) with the 
FCRAO. In order to study substruct ure with Spitzer, we have 
selected a sample of targets from the lRagan et al.l (l2006h sam- 
ple which are compact, typically 2' x 2' (or 2 x 2 pc at 
4 kpc), and opaque, providing the starkest contrast at 8 yum 
(MSX Band A) with which to examine the absorbing struc- 
ture. The selected objects also exhibit significant emission 
in transitions of CS and N2H^ that are known to trace high- 
density gas, based on their high critical densities. By se- 
lecting objects with strong emission in these lines, we en- 
sure that their densities are >10'* cm"^ and their temperatures 
are less than 20 K. Under these conditions in local clouds, 
N2H'^ is strongest when CO is depleted in the pre-stellar phase 
(iBergin & Langedll997h . hence a high N2H+/CO ratio guided 
our attempt to select the truly "starless" dark clouds in the 
IRDC sample. Our selection criteria are aimed to isolate ear- 
liest stages of star formation in local clouds and give us the 
best hope of detecting massive starless objects. The eleven 
IRDCs observed are listed in Table 1 with the distances de- 
rived in Ragan et al. (2006) using a Milky Way rotation curve 
model (iFich et al.lll989l) assuming the "near" kinematic dis- 
tance. The listed uncertainties in Table 1 arise from the +14% 
maximal deviation inherent in the rotation curve model. 



2.2. Spitzer Observations &- Data Processing 

Observations of this sample of objects were made on 2005 
May 7-9 and September 15-18 with IRAC centered on 
the coordinates listed in Table 1 . Each region was observed 

10 times with slightly offset single points in the 12s high dy- 
namic range mode. All four IRAC bands were observed over 
7'x 7' common field-of-view. MIPS observations were ob- 
tained on 2005 April 7 - 10 of the objects in this sample. Us- 
ing the "large" field size, each region was observed in 3 cycles 
for 3s at 24 yum. MIPS observations cover smaller 5.5'x 5.5' 
fields-of-view but big enough to contain the entire IRDC. Fig- 
ures [THTT] show each IRDC field in all observed wavebands. 
The absorbing structures of the IRDCs are most prominent at 
8 fim and 24 jjm. 

We used IRAC images processed by the Spitzer Science 
Center (SSC) using pipeline version S 14.0.0 to create basic 
calibrated data (BCD) images. These calibrated data were 
corrected for bright source artifacts ("banding", "pulldown", 
and "muxbleed"), cleaned of cosmic ray hits, and made 
into mosaics using Gutermuth's WCS-bas ed IRAC post- 
proce ssing and mosaicking package (see iGutermuth et al.l 
I2OO8I for further details). 

Source finding and aperture photometry were perf ormed 
using Gutermuth's Phot Vis version 1.10 ( Gutermut h et al.l 
l2008h . We used a 2.4" aperture radius and a sky annulus from 
2.4" to 6" for the IRAC photometry. The photometric zero 
points for the [3.6], [4.5], [5.8], and [8.0] bands were 22.750, 
21.995, 19.793, and 20.187 magnitudes, respectively. For the 
MIPS 24 jjm photometry, we use a 7.6" aperture with 7.6" to 
17.8" sky annuli radii and a photometric zero point of 15.646 
magnitude. All photometric zero points are calibrated for im- 
age units of DN and are corrected for the adopted apertures. 

To supplement the Spitzer photometry, we incorporate 
the source photometry from the Two-Micron All Sky Sur- 
vey (2MASS) Point Source Catalog (PSC). Source Hsts are 
matched for a final catalog by first matching the four IRAC 
band catalogs using Gutermuth's WCSphotmatch utility, en- 
forcing a 1" maximal tolerance for positive matches. Then, 
the 2MASS sources are matched with tolerance 1" to the 
mean positions from the first catalog using the same WCS- 
based utility. Finally, the MIPS 24 jjm catalog is integrated 
with matching tolerance 1.5". 

3. STELLAR CONTENT & IRDC ENVIRONMENT 

The tremendous sensitivity of Spitzer has given us the first 
abilty to characterize young stellar populations in detail. Be- 
fore the Spitzer era, IRAS led the eff'ort in identifying the 
brightest infrared point sources in the Galaxy. Only one ob- 
ject in this sample, G034.74-0.12 (Figure [TOjl has an IRAS 
point source (18526-1-0130) in the vicinity. Here, with Spitzer, 
we have identified tens of young stellar objects (YSOs) in the 
field of each IRDC. 

3.1. Young Stellar Object Identification & Classification 

With this broad spectral coverage from 2MASS to 
IR AC to MIPS, we app ly the robust critieria described 
in IGutermuth et al.l (l2008b to identify young stellar objects 
(YSOs) and classify them. Table 2 Usts the J, H, K,, 3.6, 4.5, 
5.8, 8.0 and 24 jjm photometry for all stars that met the YSO 
criteria, and we note the classification as Class I (CI), Class 

11 (CII), embedded protostars (EP), or transition disk objects 
(TD). A color-color diagram displaying these various classes 
of YSOs 
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Fig. 2. — G006.26-0.5 1 : Wavelengths as noted in Figure 1 . 
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3. — G009. 16+0.06: Wavelengths as noted in Figure 1. 
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Fig. 4. — G009.28-0.15: Wavelengths as noted in Figure 1. 
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Fig. 5. — G009.86-0.04: Wavelengths as noted in Figure 1. - Embedded 
Objects (indices 6 and 7 in Table 2 under source G009.86-0.04) are labeled. 
Source G009.86-0.04 index 6 is only detectable at 24 //m and lies right at the 
heart of the dust absorption. 
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Fig. 6. — G012.50-0.22: Wavelengths as noted in Figure 1. - Embedded 
Object (index 5 in Table 2 under source G012.50-0.22) is labeled. 
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Fig. 7. — G023.37-0.29: Wavelengths as noted in Figure 1. - Embedded 
Objects (indices 9 and 10 in Table 2 under source G023.37-0.29) are labeled. 
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Fig. 8. — G023.48-0.53: Wavelengths as noted in Figure 1. 
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Fig. 9. — G024.05-0.22: Wavelengths as noted in Figure 1. - Embedded 
Object (index 1 in Table 2 under source G24.05— 0.22) is labeled. 
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Fig. 10. — G034.74-0. 12: Wavelengths as noted in Figure 1. - Embedded 
Object (index 5 in Table 2 under source G034.74-0. 12.) is labeled. 
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Fig. 12. — IRAC four color plot for all objects in the IRDC sample for 
all objects with photometry in all four bands that had errors less than 0.2 
magnitudes. Class 1 protostars are marked with red squares, green circles 
mark the more-evolved Class 11 sources, and transition/debris disk objects 
are marked with purple circles. The deeply embedded objects identified with 
this analysis did not have sufficient detections in IRAC bands to appear on 
the color-color plots. The extinction law from . Flaherty et al. 1 2007i) indica ted 
by the black arrow, and the extinction law from ilndebetouw et alj 42003) is 
plotted as the blue arrow. 

in the entire sa mple is shown i n Figu re [T2 The extinction 



laws f rom both iFlahertv et alj (l2007h and llndebetouw et alj 
( l2005h are plotted to show the effect of five magnitudes of vi- 
sual extinction. The objects associated with these IRDCs are 
a great distance from us and in the plane of the Galaxy, so they 
naturally suffer from a great deal of extinction, reddening, and 
foreground contamination. Furthermore, the reddening law 
used in this classification scheme and the measures taken to 
extricate extragalactic contaminants may be inaccurate due to 
the great distance to IRDCs, as the criteria were originally 
designed to suit local regions. This may result in misclassi- 
fication of sources. For example, a highly reddened Class 1 
object might appear as an embedded protostar. Nonetheless, 
if these objects are indeed protostars, it is likely that they are 
associated with the IRDC. 

In Table 3, we summarize the number of each class of YSO 
in each IRDC field. We note the number of these YSOs that 
are spatially coincident with the absorbing IRDC clumps (see 
94.21 ). Only ~IQ% of the YSOs are associated with the dense 
gas. The rest appear to be a distributed population of stars 
surrounding the IRDC. This may be because any star directly 
associated with the IRDC is too heavily obscured to be de- 
tected even with the deep Spitzer observations we undertook. 
Our observatio ns are sensitive to 1 -3 Mq, IMyr-old pre-main 
sequence stars ( Baraffe et al. 1998) , or 1 Lg Class protostar 
at 4 kpc with no extinction (.Whitney et alJl2003i) . With ex- 

^ No embedde d pr otostar was detected in all four IRAC bands, so none are 
plotted in Figure [T2l 



tinction, which can reach 1 -2 magnitudes in the Spitzer bands, 
embedded YSOs up to 3-4 M© might be present, but hidden 
from our view. Another possible reason for the lack of YSOs 
detected coincident with the dense gas is that the IRDC itself 
could be in a stage prior to the onset of star formation, and the 
surrounding stars that are observed have disrupted their natal 
molecular gas. 

Table 4 lists all of the objects identified as embedded ob- 
jects that are spatially coincident with an IRDC. We list the 
flux density at each Spitzer wavelength and an estimate of the 
mid-infrared luminosity derived from integrating the spectral 
energy distribution, which is dominated by emission at 24 fim. 
In the likely event that the embedded objects are extincted, 
these mid-infrared luminosities will be underestimated. Tak- 
ing the average extinction estimations, which can be derived 
most reliably from the measurements of Class II objects. Ar- 
ranges from 1 to 3, which, if the extinction law lFlahertv et alj 
( l2007h is applied, corresponds to A24 of 0.5 to 1 .6. As a check, 
we use a second method to estimate the extinction: based on 
average values of the optical depth we measure in the IRDCs, 
we confirm that A24 ~1 is typical in these objects. Given 
the uncertain extinction properties, and the fact that a large 
portion of these embedded sources' luminosity will emerge at 
longer wavelengths not observed here, the luminosities pre- 
sented in Table 4 are lower limits. Stars with luminosities 
in this range, according to Robitaille et al. (2006), arise from 
stars ranging from 0.1 to 2 Mq, but are likely much greater 

3.2. Nebulosity at 8 and 24 

Four IRDCs in our sample (G006.26-0.51, Figure |2l 
G009. 16+0.06, Figure [2 G023. 37-0.29 Figure El 
G034.74-0.12, Figure [TOl l exhibit bright emission nebu- 
losity in the IRDC field at 8 and 24 /^m. These regions tend 
to be brightest in the thermal infrared (e.g. 24 /im) but show 
some emission at 8 /im, which suggests they are sites of high 
mass star cluster formation. To test whether the apparent 
active star formation is associated with the IRDC in question, 
or if it is in the vicinity, we correlate each instance of a 
bright emiss ion with the molec ular observations of the object 
obtained bv lRagan et alj 12006). The molecular observations 
provide velocity information which, due to Galactic rotation, 
aid in estim ating the distance to the mid-infrared emission 
dFich et al.iri989 ). This distance compared with the distance 
to the IRDC enables us to discern whether the IRDC and 
young cluster are at the same distance or one is in the 
foreground or background. 

In the case of G006.26-0.51 (Figure|2]), we detect infrared 
emission at 24 east of the IRDC. This is spatially coin- 
cident and has similar morphology to C'^O (1 -0) emission 
emitting at a characteristic velocity of 17 km s"' (Raga n et alJ 
2006), corresponding to a distance of about 3+0.5 kpc. The 
IRDC has a velocity of 23 km s which gives a distance of 
3.8 kpc, but with an un certainty of over 500 pc (see Table 1 
and lRagan et alJ (l2006l) ). Given the errors inherent in the dis- 
tance derivation from the Galactic rotation curve, we cannot 
conclusively confirm or rule out association. G009. 16+0.06 
(Figure [3]), has neither distinct velocity component evident 
in the molecular observations nor does the molecular emis- 
sion associated with the IRDC overlap with the 24 jjm emis- 
sion. Embedded clusters should be associated with molecular 
emission especially C'^O which is included in the FCRAO 
survey. Associated emission for this object likely lies out- 
side the bandpass of the FCRAO observations and is at a 
greater or lesser distance than IRDC. The 24 jum image of 
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G023.37-0.29 (Figure|7]i shows bright emission to the south 
of the IRDC and another region sHghtly south and west of the 
IRDC. This emission is not prominent in the IRAC images, 
suggesting that this is potentially an embedded star cluster 
Molecular observations show strong emission peaks in both 
CS (2-1) and N2H+ (1-0) in the vicinity of the IRAC 8 yum 
and MIPS 24 fim emission. However, there are three distinct 
velocity components evident in the observed bandpass, none 
of which is more spatially coincident with the 24 fim emis- 
sion than the others. Unfortunately, the spatial resolution of 
the FCRAO survey is insufficient for definitive correlation. 
Finally, in G34.74-0.12 (Figure fTTTi. no molecular emission 
is distinctly associated with the nebulosity; the most likely 
scenario for this object is that the associated molecular emis- 
sion lies outside the bandpass of the FCRAO observation and, 
therefore, is not associated. 

3.3. Summary of Stellar Content and IRDC Environment 

We have characterized the star formation that is possibly 
associated with the IRDCs to the extent that the Spitzer and 
millimeter data allow. The YSO population is distributed, and 
only a handful of objects identified are directly spatially as- 
sociated with the IRDC. More explicitly, in this sample, half 
(5/11) of the sample shows no clear evidence for embedded 
sources in the dense absorbing gas, and instead appear pop- 
ulated sparsely with young protostars, the photometric prop- 
erties of which are given in Table 2, and the overall IRDC 
star content is summarized in Table 3. Among those embed- 
ded objects correlated with the absorbing structure at 8 yum, 
which are summarized in Table 4, we find a marked lack of 
luminous sources (>5 L©) at these wavelengths. There may 
be significant extinction at 24 yum, in which case we would 
underestimate their luminosity. Further, even in IRDCs with 
embedded protostars, most of the cloud core mass is not as- 
sociated with an embedded source. It is our contention that 
most of the IRDC mass does not harbor significant massive 
star formation, and, hence IRDCs are in an early phase of 
cloud evolution. 

Bright emission nebulosity is evident at 8 /zm and 24 yt^m in 
four fields, presumably due to the presence of high mass stars 
or a cluster If the IRDC were associated with the nebulosity, 
it would be a strong indication that the IRDCs have massive 
star formation occurring already in the vicinity. Molecular 
data give no definitive clues that these regions are associated 



with the IRDCs. 

Most studies including this one focus primarily on the dense 
structures that comprise infrared-dark clouds, yet their con- 
nection to the surrounding environment has not yet been dis- 
cussed in the literature. While it is clear that some star for- 
mation is directly associated with the dense material, star for- 
mation is also occurring beyond the extent of the IRDC as it 
appears in absor ption. Figure [T3] shows molecular line con- 
tours from Raga n et al.l (12006) over the Spitzer 8 yum image. 
N2H+, a molecule known to trace very dense gas, corresponds 
exclusively to the dark cloud. On the other hand, C"^0 and, to 
a greater extent '^CO, show a much more extended structure, 
which suggests that the infrared-dark cloud resides within a 
greater molecular cloud complex. For all of the objects in our 
sample, the '^CO emission was present at the edge of the map 
(up to 2' away from the central position), so it is likely that the 
emission, and therefore the more diffuse cloud that it probes, 
extends beyond the mapped area. Thus, the full extent of the 
surrounding cloud is not probed by our data. 

4. TRACING MASS WITH DUST ABSORPTION AT 8 yum 

Each infrared-dark cloud features distinct absorbing struc- 
tures evident at all Spitzer wavelengths, but they are most pro- 
nounced at 8 yt^m and 24 yt^m due to strong background emis- 
sion from polycyclic aeromatic hydrocarbons (PAHs) and 
small dust grains in the respective bandpasses (Draine 2003). 
The IRDCs in this sample exhibit a range of morphologies and 
surrounding environments. Figures fTTflTlshows a morpholog- 
ical mix of filamentary dark clouds (e.g. G037.44-i-0.14, Fig- 
urefTTTi and large "round" concentrations (e.g. G006.26-0.51, 
Figure |2]i. Remarkably, these detailed structures correspond 
almost identically between the 8 yum and 24 yum bands, despite 
the fact that the source of the background radiation arises from 
separate mechanisms. At 8 jum emission from PAHs domi- 
nate on average, and at 24 //m, the bright background is due 
to the thermal emission of dust in the Galactic plane. Consid- 
ering this scenario, it is unlikely that we are mistaking random 
background fluctuations for dense, absorbing gas with the ap- 
propriate characteristics to give rise to massive star and cluster 
formation. 

4.1. Modeling the Foreground and Background 

In the Galactic plane, the 8 //m background emission varies 
on scales of a few arcminutes. To accurately estimate struc- 
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Fig. 14. — Upper left: Original IRAC 8 pm image of G024.05-0.22. Upper right: Background model using the spatial median filtering technique with a 3' 
radius. The dark cloud is virtually eliminated from the background, but still accounts for the large-scale variations. Lower left: Same as upper right panel, except 
that a 1' radius is used, which models the dark cloud as part of the background. Lower right: Same as upper right and lower left panels, except that a 5' radius is 
used, which misses the background variation and is almost a constant value. 



tures seen in absorption, we account for these variations using 
a spatia l median filte r ing tec hnique, motivated by the methods 
used in ISimon et al.l (120061) . For each pixel in the IRAC im- 
age, we compute the median value of all pixels within a vari- 
able radius and assign that value to the corresponding pixel 
in the background model. Figure [14] illustrates an example 
of several trials of this method, including models with 1', 3', 
and 5' radius of pixels included in a given pixel's median cal- 
culation. We select the size of the filter to be as small as 
possible such that the resulting map shows no absorption as 
background features. If the radius is too small, most of the in- 
cluded pixels will have low values with few representing the 
true background in the areas where absorption is concentrated 
(lower left panel of Figure[T4li. The background variations are 
also not well-represented if we select a radius too large (lower 
right panel of Figure [T4ll. Based on our analysis, the best size 
for the filter is 3'. The observed 8 //m emission is a combina- 
tion of both background and foreground contributions. 

J jestimate^^^ j l'ruedA + j I pad A (1) 

where J dA is the intensity that we measure from the 

method described above, J Ig^ dA is the true background 
intensity, which can only be observed in conjunction with 
J Ipc dA, the foreground intensity, all at 8 /zm. The relative 
importance of the foreground emission is not well-known. For 
simplicity, we assume the foreground can be approximated by 
constant fraction, x, of the emission across each field. 



J IfcdA^x j dA (2) 

One way to estimate the foreground con t ributio n has al- 
ready been demonstrated bv Johnst one et al.] (l2003h . The au- 
thors compare observations of IRDC GO 1 1.1 1-0.12 with the 
Midcourse Science Experiment (MSX) at 8 yum and the Sub- 
millimeter Common-User Bolometer Array (SCUBA) on the 
James Clerk Maxwell Telescope (JCMT) at 850 fim (see their 
Figure 3) and use the point at which the 8 /im integrated flux is 
at its lowest at high values of 850 yum flux for the foreground 
estimate . The top panel F igurefTSlshows a similar plot to Fig- 
ure 3 in iJohnstone et alJ d2003 ). except our integrated 8 /im 
flux is measured with Spitzer and presented here in units of 
MJy/sr SCUBA 850//m data for two of the IRDCs in this 
sample (G009.86-0.04 and G012.50-0.22) are available as 
part of the legacy data release (Di Francesco et al. 2008) and 
are included in this plot. Just as Johnstone et al.. (2003) point 
out, we see a clear trend: where 8 fim emission is low along 
the filament, the 850jt/m flux is at its highest. In the case 
of GO 11. 11 -0.1 2, where the SCUBA data are of the high- 
est quality, we take the minimum 8 /im flux density to be an 
estimate of the foreground contribution. Assuming this trend 
is valid for our sample of IRDCs, we use the 8 fim emission 
value measured at the dust opacity peak in each source as our 
estimation of the foreground level for that object (for the re- 
mainder of this paper, we will refer to this method as fore- 
ground estimation method "A"). Given these considerations, 
we find values for x to range between 2 and 5. Up to 20% of 
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Fig. 15.— Spitzer 8 vs. SCUBA 850pm flux for IRDC GOI 1.1 1-0.12, 
G009.86-0.04, and G012.50-0.22. The horizontal dashed line marks where 
the 8 /jm flux density reaches a minimum in GOI 1.1 1-0. 12, which is also 
indicated for the two other IRDCs with available SCUBA data. This flux 
density serves as an estimate of the foreground emission at 8 fim. The dash- 
dotted hne indicates the mean 8 fjm emission. 

this foreground contamination is likely due to scattered light 
in the detector (S.T. Megeath, private communication). We 
assume constant foreground flux at this level. As an alter- 
native foreground estimate, we also test a case in which we 
attribute half of the model flux to the background and half to 
the foreground. This is equivalent to choosing a value of x of 
1, and based on Figure[T5] is also a reasonable estimate. This 
method will be referred to as foreground estimation method 
"B." For most of the following figures and discussion, we use 
estimation method A and refer the results from method B in 
the text when applicable. 

With an estimation of the foreground contribution, the ab- 
sorption can be quantitatively linked to the optical depth of the 
cloud. The measured integrated flux, J I,„dA at any point in 
the image, including contributions from both the foreground 
and background, can then be expressed as 

J I„,dA = J I'^'^'e-'^dA + J IpcdA (3) 

where rg is the optical depth of the absorbing material. For the 
subsequent calculations, we use the average intensity, assum- 
ing uniform transmission over the IRA C channel 4 passband, 
and a verage over the extinction law ( Weingartner & Draind 
I2OOII see Section |4.2| i in this wavelength region in order to 
convert the optical depth into a column density (see discus- 
sion in next section). We note that we make no attempt to cor- 
rect for the spectral shape of the the dominant PAH emission 
feature in the 8 yum Spitzer bandpass, which we assume dom- 
inates the background radiation. In addition, clumpy material 
that may be optically thick and is not resolved by these ob- 
servations will cause us to underestimate the column density. 
These factors could introduce an uncertainty in the conver- 
sion of optical depth to column density. Still, we will show in 
Section [4.4. II that dust models compare favorably to our esti- 
mation of the dust absorption cross section, lending credence 
to our use of r as a tracer of column density. 

4.2. Identification of Structure 



Figure [16] shows a map of optical depth G024.05-0.22. 
This provides an example of the the absorbing substructure 
in one of the IRDCs in our sample. Owing to the high spa- 
tial resolution of Spitzer at 8 jim (1 pixel - 0.01 pc at 4 kpc, 
accounting for oversampling), we see substructures down to 
very small scales (~0.03 pc) in all IRDCs in our sample. 

In order to identify independent absorbing structures in the 
8 fim optic al depth map, w e employed the clumpfind algo- 
rithm (Willi ams et al.|[l994b . In the two-dimensional version, 
cl£ind2d, the algorithm calculates the location, size, and the 
peak and total flux of structures based on specified contour 
levels. We use the Spitzer PET"* to calculate the sensitivity of 
the observations, i.e. to what level the data permit us to dis- 
cern true variations from noise fluctuations. At 8 fim, the ob- 
servations are sensitive to 0.0934 MJy/sr which, on average, 
corresponds to an optical depth sensitivity (10-cr) of ~0.02. 
While the clumps take on a variety of morphologies, since 
clumpfind makes no assumptions about the clump shapes, 
we approximate the clump "size" by its effective radius, 

jnpix Apix 
reff = A —7 — 

\ ^Jos 

where n^jv is the number of pixels assigned to the clump by 
clumpfind, and Ap,., is the area subtended by a single pixel. 
The correction factor for oversampling, accounts for the 
fact that the Spitzer Space Telescope has an angular resolution 
of 2.4" at 8 yum, while the pixel scale on the IRAC chip is 1 .2", 
resulting in oversampling by a factor of 4. 

The number and size of structures identified with 
clumpfind varies depending on the number of contouring 
levels between the fixed lower threshold, which is set by the 
sensitivity of the observations, and the highest level set by the 
deepest absorption. We set the lowest contour level to lOo" 
above the average background level. In general, increasing 
the number of contour levels serves to increase the number of 
clumps found. In all cases, we reach a number of levels where 
the addition of further contouring levels results in no addi- 
tional structures. We therefore select the number of contour 
levels at which the number of clumps levels off, i.e. when the 
addition of more contour levels reveals no new clumps. We 
also remove those clumps found at the image edge or border- 
ing a star, as the background estimation is likely inaccurate 
and/or at least a portion of the clump is probably obscured by 
the star, rendering any estimation of the optical depth inaccu- 
rate. 

Using clumpfind, each IRDC broke down into tens of 
clumps, ranging in size from tens to hundreds of pixels per 
clump. The average clump size is 0.04 pc. Typically, there is 
one or two central most-massive clumps and multiple smaller 
clumps in close proximity. In some instances, clumps are 
strung along a filamentary structure, while in other cases, 
clumps are radially distributed about a highly-concentrated 
center. Figure [T7] shows an example of how the clumps are 
distributed spatially in G024.05-0.22 as clumpfind identi- 
fies them. 

With reliable identification of clumps, we next calculate in- 
dividual clump masses. As described, clumpfind gives to- 
tal optical depth measured at 8 /urn, thjoi, within the clump 
boundary, its size and position. This can be directly trans- 
formed into N(H),ot via the relationship 

^ http://ssc.spitzer.caltech.edu/tools/senspet/ 
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NiH)ro, = 



O-g fos 



(5) 



where erg is the dust absorption cross section at 8 /im. We 
derive an average value of erg over the IRAC channel 4 band- 
pass using dust models that take into account higher val- 
ues of Ry corresponding to d ense regions in the ISM. Using 
IWeingartner & Draind (1200 Ih . we use Rv — 5.5, case B val- 
ues, w hich agree with recent results from llndebetouw et al.l 
(12005). We find the value of erg to be 2.3xl0"2^cm2. 

The column density can then be used with the average 
clump size and the known distance to the IRDC, assuming 
all clumps are at approximately the IRDC distance, to find the 
clump mass. The mass of a clump is given as 



Mclwnp = 1.16mHN(H),o,Aclump 



(6) 



where nifj is the mass of hydrogen, N(H),£„ is the total col- 
umn density of hydrogen, the factor 1.16 is the correction for 
helium and Adump is the area of the clump. Table 6 gives the 
location, calculated mass and size of all the clumps identified 
with clumpfind. We also note which clumps are in the vicin- 
ity of candidate young stellar objects (Table 2) or foreground 
stars, thereby subjecting the given clump properties to greater 
uncertainty. On average (for foreground estimation method 
A), 25% of clumps border a field star, and these clumps are 
flagged and not used in the further analysis. In each infrared- 
dark cloud, we find between 3OOOM0 and IO'^Mq total mass 
in clumps, and typically -15% of that mass is found in the 
most massive clump. 

We perform the same analysis on the maps produced with 
foreground estimation method B. The foreground assumption 
in this case leads to lower optical depths across the map. 
Due to the different dynamic range in the optical depth map, 
clumpfind does not reproduce the clumps that are found with 
method A exactly. The discrepancy arises in how clumpfind 
assigns pixels in crowded regions of the optical depth map, so 
while at large the same material is counted as a clump, the ex- 
act assignment of pixels to specific clumps varies somewhat. 
On average, the clumps found in the "method B" maps tend 
to have lower masses by a factor of 2, though the sizes do not 
differ appreciably from those found with foreground estima- 
tion method A. 

4.3. Resolving Inaccuracy in Clump Mass Calculation 

The clumps identified in this fashion include a contribu- 
tion from the material in the surrounding envelope. As 
a result, a portion of the low-mass clump population may 
not be detected, and the amount of material in a given 
clump may be overestimated. To examine this effect, we 
use the gaussclumps algorithm (Stutzki & Guesten 1990) to 
identify clumps while accounting for the contribution from 
the cloud envelope. This method was designed to decom- 
pose three-dimensional molecular line observations by decon- 
voloving the data into clumps fit by Gaussians. To use the al- 
gorithm here without altering the code, we fabricated a data 
cube by essentially by mimicking a third (velocity) dimen- 
sion, thus simulating three-dimensional clu mps that were all 
center ed in velocity on a singl e central plane. iMookerjea et al.l 
(|2004') and lMotte et al.l(l2003l) have used similar techniques to 
simulate a third dimension to their dust continuum data sets. 
The gaussclumps algorithm inherently accounts for an ele- 
vated baseline level, which can be used to approximate the 
envelope. Applied to our data set, gaussclumps finds that 
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Fig. 16. — G024.05-0.22. 8 fim optical depth with contours highlighting 
the .structures. 
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Fig. 17.— G024.05-0.22. Results of the clumpfind algorithm plotted 
over Spitzer 8 fim image. Absorption identified as a "clump" is denoted by a 
number. 

15-50% of the material is in the envelope. Further discussion 
of the envelope contribution, including its eff'ect on the mass 
function, is given in 35.31 

The clumpfind and gaussclumps methods result in 
nearly one-to-one clump identification in the central region 
of the IRDC. However, because the contribution from the 
cloud envelope falls off further away from the central con- 
centration of mass in the IRDC, gaussclumps fails to find 
low-mass clumps on the outskirts of IRDCs as successfully as 
clumpfind, despite being statistically valid relative to their 
local background. We conclude that gaussclumps is not 
suitable to identify of structure in the outskirts of the IRDCs 
where the envelope is below the central level. 

Another method commonly employed in the literature to 
account for the extended structures in which dense cores re- 
side is a "wavelet subtraction" technique, which is described 
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Fig. 18. — Original optical depth image of G024.05-0.22 (left) and the wavelet subtracted image (right) of the same region. 



in lAlves et alj (120071) . To address the varying levels of back- 
ground across the optical depth map, we use the wavelet trans- 
form of the image to extract the dense cores. For one IRDC in 
our sample, G024.05-0.22, we (with the help of J. Alves, pri- 
vate communication) perform the wavelet analysis on the op- 
tical depth map. Figure [18] shows a comparison between the 
original optical depth map and the wavelet-subtracted map. 
With the removal of the "envelope" contribution in this fash- 
ion, the clumps are up to 90% less massive on average, and 
their average size decreases by 25%, or ~0.02 pc. 

Both using gaussclumps and applying wavelet subtraction 
methods to extract clumps show that the contribution of the 
cloud envelope is not yet well-constrained quantitatively. Not 
only is the cloud envelope more difficult to detect, its struc- 
ture is likely not as simple as these first order techniques have 
assumed in modeling it. As such, for the remainder of the 
paper, we will not attempt to correct the clump masses on an 
individual basis, but rather focus our attention on the clump 
population properties as a whole. In ^4.41 we employ several 
techniques to calibrate our mass estimation methods. We will 
show in ^that the effect of the envelope is systematic and 
does not skew the derived relationships, such as the slope of 
the mass function. 



4.4. Validating 8 jim absorption as a Tracer of Mass 

In previous studies, molecular clouds have been predomi- 
nantly probed with using the emission of warm dust at sub- 
millimeter wavelengths. While there are inherent uncertain- 
ties in the conversion of flux density to mass, the emission 
mechanism is well-understood. The method described above 
is a powerful way to trace mass in molecular clouds. To un- 
derstand the extent of its usefulness, here we validate dust 
absorption as a mass tracer by drawing comparisons between 
it and results using more established techniques. First, we re- 
late the dust absorption to dust emission as probes of column 
density. Second, we use observations of molecular tracers of 
dense gas not only to further cement the validity of the ab- 
sorbing structures, but also to place the IRDCs in context with 
their surroundings. Finally, we show that the sensitivity of the 
technique does not have a strong dependence on distance. 



4.4.1. Probing Column Density at Various Wavelengths 

As we discussed in ^4.11 there is an excellent correlation be- 
tween the 8 jjm and 850 /^m flux densities in IRDC GOl 1.1 1- 
0.12. Figure [TSl shows the point-to-point correlation between 
the SCUBA 850 yum flux density and Spitzer 8 /^m flux den- 
sity. This correspondence itself corroborates the use of ab- 
sorption as a dust tracer In addition, the fit to the correlation 
can confirm that the opacity ratio, /fg Mso, is consistent with 
dust behavior in high density environments. Relating the 8 jiva 
flux density 



h = fhge 



+ ffs 



(7) 



where Kg, is the 8 jum dust opacity, is the mass column 
density of emitting material, and fhg and ffg are the back- 
ground and foreground flux density estimates, respectively 
(from ^4.11 ). and the 850 fim flux density 



/850 = Bs5Q(Td = UK)Kg5m^l 



(8) 



where Bgso is the Planck function at 850 yum evaluated for a 
dust temperature of 13 K, k^sq is the dust opacity at 850 fim 
and Q is the solid angle subtended by the JCMT beam at 
850 jum, one can find a simple relation between the two by 
solving each for S(x) and equating them. The opacity ratio, 
put in terms of the flux density measurements is as follows: 



^8 
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/850 \f&-ffg 



fhg 



(9) 



From our data, we confirm this ratio is considerably lower 
(~500) in cold, high densit y environments flian in the diffuse 
interstellar dust as found bv lJohnstone et al.) (I2003h . 

We perform another consistency check between our data 
and dust models. With maps at both 8 and 24 /urn, both show- 
ing significant absorbing structure against the bright Galactic 
background (albeit at lower resolution at 24 jum), we can cal- 
culate the optical depth of at 24 /im in the same way we did 
in Section |4n The optical depth scales with the dust opacity 
by the inverse of the column density (t j oc k^/N(H)), so the 
ratio of optical depths is equal to the dust opacity ratio. We 
find that the typical ratio as measured by Spitzer in IRDCs is 
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which is comparable to 1.6, the lWeinpartner & Draind (1200 Ih 
prediction (for Ry - 5.5, ca se B) and 1-1.2 predicted by 
lOssenkopf & HenningI (1 19941) in the high-density case. We 
conclude that the dust properties we derive are consistent with 
the trends that emerge from models of dense environments 
typical of infrared-dark clouds. 

4.4.2. Molecular Line Tracers 

Molecular lines are useful probes of dense clouds, with par- 
ticular molecules being suited for specific density ranges. For 
instance, chemical models show that NaH '*' is an excellent 
tracer of dense gas in pre-stellar objects (iBergin & Langed 
Il997h . In su pport of these models, observations of l ow-mass 
dense cores (iT afalla et all 120021: iBergin et al.ll2002h demon- 
strate that N2H^ highlights regions of high central density 
(n~10'' cm"^), while CO readily freezes out onto cold grains 
(when n > 10"* cm"^), rendering it undetectable in the central 
denser regions of the cores. CO is a major destroyer of N2H^, 
and its freeze-out leads to the rapid rise in N2H^ abundance in 
cold gas. When a star is born, the CO evapo rates frorn grain s 
and N2H^ is destroyed in the proximate gas ( Lee et al]|2004l) . 
Thus, N2H^ is a preferential tracer of the densest gas that has 
not yet collapsed to form a star in low-mass pre-stellar cores. 

While N2H^ has been used extensively as a probe of the 
innermost regio ns of local cores , where densities can reach 
lO^cm"^ (e.g. Tafalla et alJl200 4'). this chemical sequence has 
not yet been observationally proven in more massi ve star 
forming regions. Noneth eless recent surveys (e.g. Sa kai et alj 
120081: iRaganet al l 120061) confirm that N2H + is prevalent in 
IRDCs, and mapping by lRagan et aP (l2006l) shows that N2H+ 
more closely follows the absorbing gas than CS or C'^O, 
which affirms that the density is sufficient for appreciable 
N2H^ emission. These single dish surveys do not have suffi- 
cient resolution to confirm the tracer's reliability on the clump 
or pre-stellar core scales in IRDCs. Interferometric obser- 
vations will be needed to validate N2H^ as a probe of the 
chemistry and dynamics of individual clumps (Ragan et al., 
in prep.). 

For one of the objects in our sample, G012.50-0.22, we had 
previous BIMA observations of N2H^ emission with 8"x4.8" 
spatial resolution. The BIMA data were reduced using the 
standard MIRIAD pipeline reduction methods (Saultetal. 
I1995I) . As in nearby clouds, such as lWalsh et alj (12004 ), the 
integrated intensity of N2H^ relates directly to the dust (mea- 
sured here in absorption) in this infrared-dark cloud. Fig- 
ure[T9]illustrates the quality of N2H^ as a tracer of dense gas, 
both in the N2H^ contours plotted over the 8 fim Spitzer im- 
age and the point-to-point correlation between the 8 //m opti- 
cal depth and the integrated intensity of N2H^. The points that 
lie above the average line, with high integrated intensities but 
low optical depth, are all in the vicinity of a foreground star in 
the 8 jjm image, which lowers our estimate for optical depth. 
In the sample, however, we have shown that the foreground 
and young stellar population is largely unassociated with the 
absorption. 

Two trends are apparent in Figure[T9l First, below t < 0.25 
there is a lack of N2H^ emission. This suggests that the ab- 
sorption may be picking up a contribution from a lower den- 
sity extended envelope that is incapable of producing signifi- 
cant N2H^ emission. This issue is discussed in greater detail 



in ^5.31 Alternatively, the interferometer may filter out ex- 
tended N2H^ emission. The second trend evident in Figure[T9l 
is that for t > 0.25, there is an excellent overall correlation, 
confirming that mid-infrared absorption in clouds at distances 
of 2 to 5 kpc is indeed tracing the column density of the dense 
gas likely dominated by pre-stellar clumps. 

In addition to directly tracing the dense gas in IRDCs, 
molecular observations can be brought to bear on critical 
questions regarding the use of absorption against the Galactic 
mid-infrared background and how best to calibrate the level 
of foreground emission. One way to approach this is to use 
the molecular emission as a tracer of the total core mass and 
compare this to the total mass estimated from 8 jim absorp- 
tion with differing assumptions regarding the contributions 
of foreground and background (see § 14.11 ). In Raga n et all 
(p006) we demonstrated that the distribution of N2H^ emis- 
sion closely matches that of the mid-infrared absorption (see 
also § lU. This is similar to the close similarity of N2H^ 
and dust continuum em ission in local pre-stellar cores (e.g. 
IBergin & Tafallair2007 '). Thus we can use the mass estimated 
from the rotational emission of N2H^ to set lim its on viable 
models of the foreground. In lRagan et al.l (l2006l) we directly 
computed a mass using an N2H^ abundance assuming local 
thermodynamic equilibrium (LTE) and using the H2 column 
density derived from the MSX 8 //m optical depth. However, 
this estimate is highly uncertain as the optical depth was de- 
rived assuming no foreground emission, and the N2H^ emis- 
sion may not be in LTE. Instead, here, we will use chemical 
theory and observations of clouds to set limits. 

N2H^ appears strong in emission in dense pre-stellar gas 
due to the freeze-out of CO, its primary destruction route. 
Detailed theoretical models of this process in g as with den- 
sities in excess of 10^ cm"^ (lAikawa et al.ll2005l) . as expected 
for IRDCs, suggest a typical abundan c e should be ~10~'° 
with respect to H2 dMaretet all |2006': 'Aikawa et al] 120051: 
Pagani et al. 200^. This value is consistent wi th that mea- 
sured in den se gas in several starless cores (iTafalla et alJ2002l: 
iMaret et ani200 6). Using this value we now have a rough test 
of our foreground and background estimates. For example, in 
G024.05-0.22 we find a total mass of 4100 Mq (foreground 
estimation method A). Using the data in Ragan et al. (2006), 
we find that the total mass traced by N2H^ is 4400 Mq, pro- 
viding support for our assumptions. Figure |20] shows the re- 
lationship between the total clump mass derived from absorp- 
tion and the total mass derived from our low-resolution maps 
of N2H^ for the eight IRDCs in our sample that were detected 
in N2H^. In general, there is good agreement. We plot a 30% 
systematic error in the total clump masses (abscissa) and a 
factor of 5 in for the total N2H^ mass estimate (ordinate). In 
the cases where the estimates differ, the N2H^ mass estimate 
tends to be greater than the total mass derived from the dust 
absorption clumps. This discrepancy likely arises in large part 
from an under-estimation of N2H^ abundance and/or non- 
LTE conditions. All the same, the consistency of the mass 
estimates, together with the morphological correspondence, 
reaffirms that the we are probing the dense clumps in IRDCs 
and that our mass probe is reasonably calibrated. 

We find no discernible difference between methods A and 
B of foreground estimation. However, we note that both are 
substantially better than assuming no foreground contribution. 
We therefore believe that method A is an appropriate estimate 
of the foreground contribution (see ^4.11 ). 

4.4.3. Effects of Distance on Sensitivity 
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Fig. 19. — Left: Contours of integrated intensity of N2H^ plotted over the IRAC 8 /im image of IRDC G012.50-0.22. Right: Point-to-point correlation of the 
N2H^ integrated intensity and the 8 fjm optical depth. Points with high integrated intensity but low optical depth correspond to stars, whose presence leads to the 
underestimation of optical depth in the vicinity. 
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Fig. 20. — Comparison of the total mass derived from NaH^ maps from 
IRagan et alj ^006) and total clump mass as derived from dust absorption at 
8 ^m, where the black diamonds represent the mass using foreground esti- 
mation method A and the grey squares show the masses derived using fore- 
ground estimation method B (see i|4.U . Thi'ee of the IRDCs in the sample 
did not have adequate NtH^ detections. Error bars for 30% systematic er- 
rors in the mass ai'e plotted for the clump mass estimates, and a factor of 5 
uncertainty is plotted for the N2H+ mass estimates. The dashed line shows a 
one-to-one correspondence for reference. 

Infrared-dark clouds are much more distant than the lo- 
cal, well-studied clouds such as Taurus or p Ophiuchus. As 
such, a clear concern is that the distance to IRDCs may pre- 
clude a well-defined census of the clump population. The 
most likely way in which the our survey is incomplete is the 



under-representation of low-mass objects due to their rela- 
tively small size, blending of clumps along the line of sight, or 
insensitivity to their absorption against the background. One 
observable consequence of this effect, assuming IRDCs are a 
structurally homogeneous class of objects, might be that more 
distant IRDCs should exhibit a greater number of massive 
clumps at the expense of the combination of multiple smaller 
clumps. Another possible effect is the greater the distance 
to the IRDC, the less sensitive we become to small clumps, 
and clumps should appear to blend together (i.e. neighboring 
clumps will appear as one giant clump). Due to this effect, we 
expect that the most massive clumps of the population will be 
over-represented. As a test, we examine the distribution of 
masses and sizes of clumps as a function of IRDC distance, 
which is shown in Figure [2T| This sample, with IRDCs rang- 
ing in distance from 2.4 to 4.9 kpc away, does not show a 
strong trend of this nature. We show the detection limit for 
clumps to illustrate the very good sensitivity of this technique 
and that while it does impose a lower boundary on clump de- 
tectability, most clumps are not close to this value. We found 
no strong dependence of clump mass or size on the distance 
to the IRDC and conclude that blending of clumps does not 
have a great effect on the mass sensitivity. 

Typical lo w-mass star forming cores range in size from 0.03 
to 0. 1 pc (Bergin & Tafallal2007 ). If one were to observe such 
objects 4 kpc, they would only subtend a few arcseconds. For 
example, if LI 544, a prototypical pre-stellar core, resided at 
the typical distance to the IRDCs in the sample, it would show 
sufficient abs orption (based on rep orted column density mea- 
surements by lBacmann et al."2000') against the Galactic back- 
ground, but according to 'Williams et al. (2006), would sub- 
tend 3" in diameter at our fiducial 4 kpc distance, which is 
very close to our detection limit. In addition, very low mass 
clumps could blended into any extended low-density material 
that is included in our absorption measurement. These effects 
should limit our sensitivity to the very low-mass end of our 
clump mass function. 

To first order, we have shown distance is not a major factor 
because the high-resolution offered by Spitzer improves our 
sensitivity to small structures. However, infrared-dark clouds 
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Fig. 21. — Top: The range in clump mass as a function of distance. The 
median clump mass for each IRDC in the sample is indicated with a diamond. 
Bottom: The range in clump size as a function of distance. The median clump 
size for each IRDC in the sample is indicated with a diamond. The resolution 
limit is plotted as a solid line, and it shows the boundary at which clumpfind 
defines a "clump" for an object at the distance of the indicated host IRDC. 

are forming star clusters and by nature are highly structured 
and clustered. As such, we can not rule out significant line- 
of-sight structure. Since independent clumps along the line- 
of-sight might have distinct characteristic velocities, the addi- 
tion of kinematical information from high-resolution molecu- 
lar data (Ragan et al., in prep.) will help the disentanglement. 

5. MASS FUNCTION 

A primary goal of this study is to explore the mass function 
of clumps in infrared-dark clouds and compare it to that of 
massive star formation regions, local star formation regions, 
and the stellar IMF. We note that there is some ambiguity in 
the literature about the "clump" versus the "core" mass func- 
tions. In the following description, a "core" mass function 
refers to the mass spectrum objects with masses in the "core" 
regime (10"' -10' Mq, 10^^-10"' pc), and a "clump" mass 
function for objects in the "clum p" regime (10~'-10' Mq , 
10^^-10 ' pc), as summarized in iBergin & Tafallal (l2007h . 
Here we present the infrared-dark cloud clump mass func- 
tion. We describe the relevance of this result in the context of 
Galactic star formation and discuss several methods we use to 
test its validity. 

5.1. Mass Function in Context 

A fundamental property of the star formation process is the 
mass spectrum of stars, and, more recently, the mass func- 
tion of pre-stellar objects. The mass spectrum in either case 
is most typically characterized by a power law, taking the 
form dN/dM oc M^", known as the differential mass func- 
tion (DMF). In other contexts, the mass function can be de- 
scribed as a function of the logarithm of mass, which is con- 
ventionally presented as dN/d(log m) oc , in which case 
r - -{a - 1). In the results that follow, we present the slope 
of the clump mass function in terms of a. 

A commonly used method for studying mass functions of 
pre-stellar cores is observation of dust thermal continuum 
emission in nearby star-forming clouds. Cold dust emis- 



sion is optically thin at millimeter and sub-millimeter wave- 
lengths, and can therefore be used as a direc t tracer of mass. 
A number of surveys o f local clouds (e.g. iJohnstone et alj 
2000; Mott e et anil998h have been performed with single- 
dish telescopes covering large regions in an effort to get 
a complete picture of the mass distribution of low-mass 
clouds. This i s an extremely powerful technique, but as 
iGoodman et al.l (12008 ) demonstrate, this technique suffers 
from some Umitations, chief among them poor spatial resolu- 
tion (in sin gle-dish studies), req uired knowledge of dust tem- 
peratures ( P avlvuchenkov et al.ii200 7). and the insensitivity to 
diffuse extended structures. 

Another technique that has been empl oyed to map dust 
employs near-infrared extinction mapping (lAlves et al.ll2"007t 
Lomba rdi et al.l 120061) , which is a way of measuring Ay due 
to da rk clouds by probin g the color excesses of background 
stars jLombardi & Alvesl i2001). This method is restricted to 
nearby regions of the Galaxy because of sensitivity limita- 
tions and the intervention of foreground stars, both of which 
worsen with greater distance. Also, the dynamic range of Ay 
in such studies is limited to ~l-60 (Lombardi & Alves 2001), 
while our technique probes from Ay of a few to ~ 100. 

The dust-probing methods mentioned above, both thermal 
emission from the grains and extinction measures using back- 
ground stars, often find a core mass function (CMP) that is 
similar in shape to th e stella r initial mass function (IMF), as 
describe d by Salpeter (Il955h . where a = 2.35 (T = -1.35), or 
iKroupal (12001.) . This potentially suggests a one-to-one map- 
ping between the CMP an d IMF, perhaps scaled by a constant 
"efficiency" factor (e.g. Alves et ani2007 ). Also, both tech- 
niques are difficult to apply to regions such as infrared-dark 
clouds due to their much greater distance. As we show in 
94.21 absorbing structure exists below the spatial resolution 
limit of single-dish surveys. Sensitivity limitations and fore- 
ground contamination preclude use of extinction mapping to 
probe IRDCs. 

Structural analysis using emission from CO isotopologues 
find a somewhat diff erent character to the distribution of mass 
in molecular clouds. iKramer et al.l (1199^ determined that the 
clump mass function in molecular clouds follows a power law 
with a between 1.4 and 1.8 (-0.8 < F < -0.4). This is sig- 
nificantly shallower than the Salpeter-like slope for clumps 
found in works using dust as a mass probe. This disagree- 
ment may be due to an erroneous assumption about one tech- 
nique or the other, or it may be that the techniques are find- 
ing information about how the fragmentation process takes 
place from large scale, probed by CO, to small scales, probed 
by dus t. Another possible e xplanation is that most of the ob- 
jects in lKramer et al.l ( Il998h are massive star forming regions, 
and star formation in these regions may be intrinsically differ- 
ent than tyical regions studied in the local neighborhood (e.g. 
Taurus, Serpens). 

Sub-millimeter observations of more dista nt, massive star- 
formation regions have been undertaken (e.g.lReid & WilsonI 
l20Q6t iLi et alJl2007t iMookerjea et"an 120041: iRathborne et alJ 
2006) with a mixture of results r egarding the mass func- 
tion shape. Rathbor ne et al.l ([2006'), for example, performed 
IRAM observations of a large sample of infrared-dark clouds. 
Each cloud in that sample is comprised of anywhere from 2 
to 18 cores with masses ranging from 8 to 2OOOM0. They 
find a Salpeter-like {a ~2.35) mass function for IRDC cores. 
However, our Spitzer observations reveal sign ificant s t ructur e 
below the spatial resolution scales of Rathb orne et al.l (l2006l) . 
As we will show (see Section |5]), the mass function within a 
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Fig. 22. — Differential mass function of ensemble IRDC sample. Black 
filled circles indicate results of the clumpfind technique, and the green 
open triangles denote the results of the gaussclumps clump-finding method. 
The fits are broken power laws. On the high-mass end, the slope of the 
gaussclumps method mass function (a = 1.15 ± 0.04) is shallower than 
the slope of hte clumpfind mass function (a = 1 .76 ± 0.05). 

fragmenting IRDC is shallower than Salpeter and closer to the 
mass function derived from CO emission. 

Given the strong evidence for fragmentation, it is clear that 
IRDCs are the precursors to massive clusters. We then nat- 
urally draw comparisons between the characteristics of frag- 
menting IRDCs and the nearest region forming massive stars, 
Orion. At ~500 pc, it is possible to resolve what are likely to 
be pre-stellar objects in Orion individually with current obser- 
vational capabilities. With the high-resolution of our study, 
we can examine star formation regions (IRDCs) at a simi- 
lar level of detail as single-dish telescopes can survey Orion. 
For example, we detect structures on th e same size sc ale 
(~0.03 pc) as the quiescent cores found by iLi et al] (|2007|) in 
the Orion Molecular Cloud, however the most massive core in 
their study is ~50 Mq. These cores account for only a small 
fraction of the total mass in Orion. 

5.2. Results: Differential Mass Function 

We use the IRDC clump masses calculated in 34.21 (using 
clumpfind and foreground estimation method A) to con- 
struct an ensemble mass function in Figure |22] The mass 
function that results from using foreground estimation method 
B is shifted to lower masses by a factor of 2, but the shape is 
identical. Because IRDCs appear to be in a roughly uniform 
evolutionary state over the sample (i.e. they are all likely as- 
sociated with the Molecular Ring, and they possess similar 
densities and temperatures), we merge all the clumps listed in 
Table 6 as ensemble and present a single mass function for all 
the objects at a range of distances. This assumes that the char- 
acter of the mass function is independent of the distance to a 
given IRDC. Recall that we see no evidence (see Figure ISTT l 
for the mass distributions to vary significantly with distance. 

For the calculation of the errors in the DMF we have sep- 
arately accounted for the error in the mass calculation and 
the counting s t atistic s. We used a method motivated by 
iReid & WilsonI ( |2005|) to calculate the mass error We have 



assumed that the clump mass error is dominated by the sys- 
tematic uncertainty of 30% in the optical depth to mass cor- 
rection. For each clump we have randomly sampled a Guas- 
sian probability function within the Icr envelope defined by 
the percentage error. With these new clump masses we have 
re-determined the differential mass function. This process is 
repeated lO'* times, and the standard deviation of the DMF in- 
duced by the error in the mass is calculated from the original 
DMF. This error is added in quadrature to the error introduced 
by counting statistics. The provided errors are Icr, with the 
caveat that the value assumed for the systematic uncertainty 
is open to debate. As a result, when there are large numbers 
in a given mass bin, the error is dominated by the mass uncer- 
tainty. Conversely, when there are few objects in a mass bin, 
the error is dominated by counting. 

The IRDC clump mass function for this sample spans 
nearly four orders of magnitude in mass. We fit the mass 
function with a broken power law weighted by the uncer- 
tainties. At masses greater than ~40Mq, the mass func- 
tion is fit with a power law of slope Qf=1.76+0.05. Below 
~4QMq, the slope becomes much shallower, a=0.52+0.04. 
We also include in Figure |22] the mass function of clumps 
found with the gaussclumps algorithm, with errors calcu- 
lated in the identical fashion. Performing fits in the equiv- 
alent mass regimes results in a shallower slope for masses 
greater than 40Mq (a'=1.15±0.04), while the behavior at low 
masses is similar. As discussed in 34.21 the clumps found with 
clumpfind and gaussclumps are in good agreement in the 
central region of each IRDC, but tend to disagree on the out- 
skirts. This is a consequence of the failure of gaussclumps 
to model the varying background. Examination of the im- 
ages reveals that the contribution of the diffuse material varies 
across the image, thereby setting the background level too 
high for outer clumps (where the envelope contributes less) 
to be detected. In fact, these clumps appear to be preferen- 
tially in the 30 to SOOMq range, and a mass function con- 
structed with the gaussclumps result is significantly shal- 
lower than derived with clumpfind (see Figurel22li. We con- 
clude that gaussclumps is not suitable to identify structure 
away from the central region of the IRDC where the envelope 
level is below the central level. This is further supported by 
the wavelet analysis which is capable of accounting for a vari- 
able envelope contribution. It is worth noting that for the one 
IRDC for which we have the wavelet analysis, that the slope 
of the derived mass function shows little appreciable change 
and agrees with the clumpfind result. 

To put the mass function into context with known Galac- 
tic star formation, we plot the clump mass function of all 
clumps in our sample in Figure |23] along with the core/clump 
mass function of a number of other studies probing vari- 
ous mass ranges. We select four studies, each probing mas- 
sive star forming regions at different wavel engths an d res- 
olutions includi ng quiescent cores in Orion (iLi et alj|2007,) , 
clumps in MI7 (iReid & Wilsonl l2006h. clumps in RCW 106 
(iMookerjea et al. *2004'). and clumps in massive star formation 
region NGC 6334 (Munoz et al. 2007). In their papers, each 
author presents the mass function in a different way, making it 
difficult to compare the results directly to one another Here, 
we recompute the mass function for the published masses in 
each work uniformly (including the treatment of errors, see 
above). Each of the mass functions is fit with a power law. 
Figure |23] highlights the uniqueness of our study in that it 
spans over a much larger range in masses than any other study 
to date. 
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At the high-mass end, th e m ass function agrees w ell with 
the iMookerjea et alJ (l2004 and lMunoz et alJ (l2007h studies, 
which probed to lower mass limits of 30Mq and 4Mq, respec- 
tively. The fall-off from the steep slope at the high mass end to 
a shallower slope at the low mass end immediately suggests 
that completeness, enhanced contribution from the envelope 
and/or clump blending become an issue. Ho wever, the slope 
at t he low mass end compares favorably with lLi et alj (l2007h 
and iReid & WilsonI (2006) which probe mass ranges 0.1 to 
46Mq and 0.3 to 2OOM0, respectively. In addition to the gen- 
eral DMF shape at both the high mass and low mass end, the 
"break" in the mass function falls in the lOM© to 5QMq range 
for the ensemble of studies, including ours. If this is a real 
feature of the evolving mass spectrum, this can shed some 
light on the progression of the fragmentation process from 
large, massive objects to the numerous low-mass objects like 
we see in the local neighborhood. The characteristic "break" 
mass can also be a superficial artifact of differences in bin- 
ning, mass determination technique, and observational sensi- 
tivity. Our study is the only one that spans both mass regimes, 
and further such work is needed to explore the authenticity of 
this feature. However, in ^we speculate that this may be an 
intrinsic feature. 

It is possible that the slope of the IRDC clump mass func- 
tion might be an artifact of a limitation in our technique. With 
the great distances to these clouds, one would expect the ef- 
fect of clump blending to play a role in the shape of their mass 
spectrum. We have shown in § 14.4.31 that distance does not 
dramatically hinder the detection of small clumps. Our study 
samples infrared-dark clouds from 2.4 kpc to 4.9 kpc, and we 
find that the number of clumps does not decrease with greater 
distance, nor does the median mass tend to be be significantly 
greater with distance. Furthermore, with the present analy- 
sis, we see no evidence that including clumps from IRDCs at 
various distances affects the shape of the mass function. 

From past studies of local clouds there has been a dispar- 
ity between mass functi on slope derived with dust emission 
and C O (e.g. compare iJohnstone et al.ll2001l : iKramer et alj 
Il998l) . Our result suggests that massive star forming regions 
have mass functions with slope in good agreement with CO 
isotopologues, e.g. a-l.E. This is crucial because CO ob- 
servations contain velocity information, which allow for the 
clumps to be decomposed along the line-of-sight. Still, the 
authors find a shallow slope in agreement with ours. We con- 
clude that clump blending, while unavoidable to some extent, 
does not skew the shape of the mass function as de rived from 
dust e mission or absorption. A close look at Kram er et alj 
(119981) results finds that the majority of objects studied are 
massive star formation regions. Given the general agreement 
of the clump mass function of this sample of IRDCs with 
other studies of massive star formation regions, we believe 
this result represents the true character of these objects, not 
an artifact of the observing technique. 

Several studies of pre-stellar cores in the local neighbor- 
hood show a mass distribution that mimics the shape of the 
stellar IMF. That the slope of the mass function in infrared- 
dark clouds is considerably shallower than the stellar IMF 
should not be surprising. The masses we estimate for these 
clumps are unlikely to give rise to single stars. Instead, the 
clumps themselves must fragment further and eventually form 
a star cluster, likely containing multiple massive stars. Unlike 
Orion A, for example, which contains ~10'' Mq distributed 
over a 380 square par sec (6.2 square degrees at 450 pc) re- 
gion (ICarpentedl2000l) . in IRDCs, a similar amount of mass 
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Fig. 23. — Differential mass function of this IRDC sample (black filled cir- 
cles) fit with a single power-law for Mc/„„,p > 30Mq (a = 1.76 ± 0.05) com- 
pared with various star formation regions in the high mass regime and their 
respective single power-law fit slopes. At the high mass end, our fit agrees 
well with that of other studies: Open purple diamonds from Muiioz et^ 
t2007) (a = 1.64±0.06); Open green inverted triangles from Mookerj ea etaU 
^004) (a = 1.59 ± 0.10). At the low mass end, we fit a second power law 
for the bins with Mc;„„,p < 30Mo (a = 0.52 ± 0.04), which agrees well with 
other studies in this mass regime: Open blue diamonds from Reid & Wilso^ 
(2006) (a = 0.80 ± 0.07); Open red circles from quiescent Orion cores from 
Li et al. l2007)(a = 0.82 ± 0.09). Note that only this study spans the entire 
range of masses, so the reality of the apparent break at ~30Mo is in question. 

is concentrated in clumps extending only a 1.5 square parsec 
area. Therefore, we posit that IRDCs are not distant analogues 
to Orion, but more compact complexes capable of star forma- 
tion on a more massive scale. 

Given the high masses estimated for infrared-dark clouds, 
yet the lack evidence for the massive stars they must form 
is perhaps indicating that we see them necessarily because 
we are capturing them just before the onset of star formation. 
Such a selection effect would mean that we preferentially ob- 
serve these dark objects because massive stars have yet to dis- 
rupt their natal cloud drastically in the process of protostar 
formation. 

5.3. The Contribution from the IRDC Envelope 

Like nearby clouds, infrared-dark clouds are structured hi- 
erarchically, consisting of dense condensations embedded in 
a more diffuse envelope. Here we present various attempts 
to estimate the fraction of the total cloud mass resides in 
dense clumps compared to the extended clouds. First, we use 
archival '^'CO data to probe the diffuse gas and use it to esti- 
mate the envelope mass. To further explore the contribution of 
the envelope, we demonstrate that a wavelet analysis, a tech- 
nique designed to remove extended structures from emission 
maps, gives a similar relationship between envelope and dense 
clump mass. Alternatively, applying the gaussclumps algo- 
rithm to the data provides an average threshold that describes 
the diffuse structure. 

We use '^CO (1-0) molecular lin e data from the Galac- 
tic Ring Survey (iJackson et al.ll200"6 ) in the area covered by 
our Spitzer observations of G024.05-0.22 to probe the diffuse 
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Fig. 24. — Left: The mass-radius relationship for clumpfind clumps (foreground method A) in the entire sample of IRDCs (gray), with the clumps found only 
in G024.05-0.22 highlighted in black, and the clumps found in the wavelet subtracted image (red). The solid line denotes the critical Bonnor-Ebert mass-radius 
relation for T„„j„,o/=15 K. The dashed line is the M oc R^ - from the IKramer et al.l (1996) CO multi-line study of Orion. The dash-dotted line is taken from 
IWilliamsetal] 11994). which finds M cc R^'. Right: The mass-radius relationship for IRDC clumps, including a comparison to all the studies of massive star 
forming regions included in Figure |23] 



material in the field. The '^CO emission is widespread, cover- 
ing the entire area in the IRAC field, thus we are not probing 
the entire cloud. Assuming local thermodynamic equilibrium 
(LTE) at a temper ature of 15 K and a '^C O abundance relative 
to H2 of 4 X 10"^ ([Goldsmith et alj|2008h . we find that the the 
clump mass is ~2Q% of the total cloud mass. This demon- 
strates that IRDCs are the densest regions of much larger 
molecular cloud complexes; however, the fraction of mass 
that we estimate the clumps comprise relative to the cloud 
is an upper limit because the full extent of the cloud is not 
probed with these data. 

In 34.31 we discuss two ways in which we account 
for the envelope in the clump-finding process. First, the 
gaussclumps algorithm is an alternative method of identi- 
fying clumps, and in 35.2l we examine the effect this method 
has on the clump mass function. The algorithm is insensi- 
tive to clumps on the outskirts of the IRDC, thereby flattening 
the mass function. While gaussclumps may oversimplify 
the structure of the envelope for the purposes of identifying 
clumps, it does provide a envelope threshold, above which 
optical depth peaks fit as clumps and below which emission 
is subtracted. This threshold approximates the level of the 
envelope, and as a result, gaussclumps finds 15-50% of the 
optical depth level is from the diffuse envelope. The wavelet 
subtraction technique results in clumps that are on average 
90% less massive and smaller in size by 25% (~0.02 pc) than 
those extracted from the unaltered map. 

These analyses of the IRDC envelope show us that our tech- 
nique is only sampling 20-40% of the clouds total mass and, 
at the same time, the clump masses themselves include a con- 
tribution from the surrounding envelope. Because of these 
factors, the different methods for isolating "clumps" have 
varying levels of success. For example, using gaussclumps 



equips us to parametrically remove the envelope component 
to the clump, but due to the underlying assumption of the 
baseline level, it misses many clumps that clumpfind iden- 
tifies successfully. The mass function that results from us- 
ing the gaussclumps method is shallower than that from 
clumpfind, as gaussclumps fails to find clumps on the 
periphery of the dominant (often central) concentration of 
clumps, where the envelope level is lower 

While both the clumpfind and gaussclumps methods 
have their drawbacks, it is clear that IRDCs have significant 
structure on a large range of scales. The relatively shallow 
mass function for IRDC clumps and other massive star form- 
ing regions shows that there is a great deal of mass in large 
objects, and future work is needed to understand the detailed 
relationship between the dense clumps and their surroundings. 

6. MASS-RADIUS RELATION 

Next we investigate the relationship between the mass and 
size of the clumps found in IRDCs, which informs us of the 
overall stability of the clump structures. Figure |24] shows the 
mass-radius relationship of the clumpfind-identified clumps, 
highlighting the results for G024.05-0.22 and the wavelet- 
subtracted case. Indeed, the clumps extracted from the 
wavelet-subtracted map are shifted down in mass by 90% 
and down in size by 25%, but the relationship between the 
quantities does not change. We plot the relation of simple 
self-gravitating Bonnor-Ebert spheres (M (7?) - lARcp-jG, 
where a is the sound speed and set to 0.2 km s solid line) 
and also the mass-radius relationship ob served in a m ulti-line 
CO survey of Orion (M oc R^^, Ki-amer 'etaLlll996i (dashed 
line)). For comparison. Figure |24] also shows these proper- 
ties from the other studies of massive star formation regions. 
We note that the spatial resolution of the comparison studies 



26 



Ragan et al. 



is larger t han the resolu tion of this study. The relat ionship 
for Or ion (iLi et aUlIOOTh Ml 7 (iReid & WilsonllliOfih. NGC 
6334 jMunoz et all I2007h and RCW 106 jMookeriea et alj 
12004 all agree with the iKramer et all ( Il996h relationship, 
which is con siste nt with the mass function agreement to CO 
studies (see ^5. 2b . 

The IRDC clumps are clearly gravitationally unstable, 
showing higher densities than their local Bonnor-Ebert sphere 
counterparts. The relationship for clumps in IRDC s shows a 
steeper trend, one closer to the Williams et al.l ( 1199 4) relation- 
ship, M oc 7?^-^. Also, dust extinction at 8 /im has greater sen- 
sitivity to high-densities than CO, which is known to freeze 
out at extreme densities. Hence, while the IRDC clumps are 
clearly Jeans unstable, the slope of the relation may be simply 
a reflection of the different mass probe used here. 

7. DISCUSSION & CONCLUSION 

The Spitzer Space Telescope affords us the ability to probe a 
spatial regime of massive clouds in the Galactic Ring at com- 
parable resolution as has been applied to the numerous studies 
of local, low-mass star formation. In this way, we can extend 
the frontier of detailed star formation studies to include re- 
gions the likes of which are not available in the solar neighbor- 
hood. This study demonstrates a powerful method for charac- 
terizing infrared-dark clouds, the precursors to massive stars 
and star clusters. These objects provide a unique look at the 
initial conditions of star formation in the Galactic Ring, the 
dominant mode of star formation in the Galaxy. 

We present new Spitzer IRAC and MIPS 24 yum photometric 
measurements supplemented with 2MASS J, H, K, photom- 
etry of the distributed young stellar population observed in 
the Spitzer fields. Rigid color criteria are applied to identify 
candidate young stellar objects that are potentially associated 
with the infrared-dark clouds. In all, 308 young stellar objects 
were identified (see Table 2), seven of which are classified as 
embedded protostars. For those objects, we set lower limits 
on the infrared luminosities. One IRDC has an IRAS source 
in the field, which is the best candidate for an associated mas- 
sive star. Otherwise, our observations provide no evidence 
for massive star formation in IRDCs, though sensitivity lim- 
itations do not rule out the presence of low mass stars and 
heavily extincted stars. Nebulosity at 8 and 24 yum was de- 
tected in four of the fields, but when these regions were corre- 
lated with molecular data, they do not appear to be associated 
with the IRDCs. On average, 25% of clumps are in the vicin- 
ity of stars and ~10% are in near YSOs, which are the most 
likely sources to be associated with the infrared-dark cloud. 
Since most of the mass is not associated with any indicator of 
star formation. This leads us to conclude that IRDCs are at 
in earlier stage than, say, the nearest example of massive star 
formation, the Orion Nebula, and these results are powerful 
clues to the initial conditions of star cluster formation. 

We detail our method of probing mass in IRDCs using dust 
absorption as a direct tracer of column density. We perform 
the analysis using two different assumptions (methods A and 
B) for the foreground contribution to the 8 jum flux. The 
IRDC envelope contribution to the To validate our method 
in the context of others, we compare and find good agree- 
ment between the 8 yum absorption and other tracers of dust, 
such as sub-millimeter emission from dust grains measured 
with SCUBA and N2H^ molecular line emission measured 
with FCRAO and BIMA. We show that distance does not play 
a role in the effectiveness of the technique. The high reso- 
lution Spitzer observations allows us to probe the absorbing 



structures in infrared-dark clouds at sub-parsec spatial scales. 
We apply the clumpfind algorithm to identify independent 
absorbing structures and use the output to derive the mass 
and size of the clumps. Tens of clumps are detected in each 
IRDC, ranging in mass from 0.5 to a few x 10^ Mq with 
sizes from 0.02 to 0.3 pc in diameter. We also apply the 
gaussclumps algorithm to identify clumps. The structures 
in the central region of the IRDC correspond almost perfectly 
to the clumpfind result, but gaussclumps misses clumps on 
the outskirts because it fails to account for a spatially variable 
background level. 

The existence of substructure - from 10^ Mq clumps down 
to 0.5 Mq "cores" - indicates that IRDCs are undergoing frag- 
mentation and will ultimately form star clusters. The typi- 
cal densities (n > 10^ cm"-') and temperatures (T < 20 K) 
of IRDCs are consistent with massive star forming regions, 
but they lack the stellar content seen in more active massive 
star formation regions, such as the Orion molecular cloud or 
W49, for example. The mass available in the most massive 
clumps, however, leads us to conclude that IRDCs will even- 
tually form multiple massive stars. 

The IRDC clump mass function, with slope a - 1 .76 + 0.05 
for masses greater than ~4OM0, agrees with the mass func- 
tion we calculate based on data from other studies of mas- 
sive objects. The mass function for both IRDCs and these 
massive clump distributions is shallower than the Salpeter- 
like core mass function reported in local regions. In fact, 
the IRDC clump mass function is more consistent with that 
found when probing molecular cloud structure using CO line 
emission {a - 1.6 - 1.8 ), again supporting the assertion that 
these objects are at an earlier phase of fragmentation. At the 
low-mass end (M < 4OM0), we find a much shallower slope, 
a - 0.52 + 0.04, which is somewhat flatter than other stud- 
ies that cover the same range in masses. This could be due 
in part to incomplete sampling of the fields. Alternatively, the 
apparent flattening of the clumps mass function around 40 
could indicate a transition between objects that will generate 
clustered star formatior i and those that give ris e to more dis- 
tributed star formation d Adams & Mversir2001h . 

IRDC clumps are generally not in thermodynamic equilib- 
rium, but rather are undergoing turbulent fragmentation. The 
mass spectrum is consistent with the predic tions of gravo- 
turbulent fragmentation of molecular clouds (lKlessenll200L) . 
The dynamic Molecular Ring environment could naturally 
be conducive for producing concentrated cluster-forming re- 
gions. 

Just as in all surveys of IRDCs to date, this study is subject 
to the blending of clumps, which could alter the shape of the 
mass function to over-represent the most massive clumps at 
the expense of clumps of all masses and sizes. To the extent 
that this sample allows, we find that this does not drastically 
effect the shape of the mass function. Other studies of cloud 
fragmentation that have the advantage of a third dimension of 
information also fin d a shallower clump mass function slope 
(iKramer et al.lll998h . We therefore conclude that this result 
is a true reflection of the structure in IRDCs and nature of 
massive star formation. 

Infrared-dark clouds are already well-established candi- 
dates for the precursors to stellar clusters and exhibit sig- 
nificant structures down to 0.02 pc scales. The properties 
of IRDCs provide powerful constraints on the initial condi- 
tions of massive and clustered star formation. We suggest 
that the mass function is an evolving entity, with infrared- 
dark clouds marking one of the earliest stages of cluster for- 
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mation. The mass distribution is top-heavy, with most of 
the mass in the largest structures. As the massive clumps 
fragment further, the mass function will evolve and become 
steeper The clumps will ultimately fragment to the stellar 
scale and then take on the Salpeter core mass function that 
has been observed so prevalently in local clouds. For exam- 
ple, following the (mostly) starless IRDC phase of cluster evo- 
lution, the mass spectrum will evolve into its steeper form, 
ahgning with the rn ass function of local embedded clusters 
(lLada& Ladall2003h or star clusters in the Large and Small 
Magellanic Clouds dHunter et al.ll2003l) . with a slope a ~ 2. 
As fragmentation proceeds on smaller scales, the mass func- 



tion would take on yet a steeper char acter observed in core 
mas s functions (e. g. lAlves et al.ll20()7h and, ultimately, stars 
(e.g. lKroupall200l . 
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TABLE 1 
Spitzer targets 



IRDC 






distsncc 




(J2000) 


(J2000) 


(kpc) 


G05. 85-0.23 
O06.26-0.51 
G09. 16+0.06 
G09.28-0.15 


17:59:53 
18:01:50 
18:05:50 
18:06:54 


-24:00:10 
—23:47:1 1 
-20:59:12 
-20:58:51 


•3 14+0.66 

O Q+0.61 

-^■°-0.69 

4,48+0.54 


G09.86-0.04 


18:07:40 


-20:25:25 


2 36+^-^^ 


G12.50-0.22 


18:13:45 


-18:11:53 


T <:c+0.67 


G23.37-0.29 


18:34:51 


-08:38:58 


4 70+0-* 
-Q.S8 


G23.48-0.53 


18:35:57 


-08:39:46 


4 10+"-^^ 
-Q.9Q 


G24.05-0.22 
G34.74-0.12 


18:35:52 
18:55:14 


-08:00:38 
+01:33:42 


4.86_i.45 


G37.44+0.14 


18:59:08 


+04:03:31 


9 sQ+1.47 
^■-^^-1.34 



" Dista nces calculated from Galactic rotation curve, as presented in lRagan et alj 



TABLE 2 

SprrzER-roENTiFiED YSOs: 2MASS, IRAC, and MIPS photometry 



1 ix ! JK. uitiex 




yjzynAf} 


T 
J 


n 




F'? fin 




p.ej 








ciassitication 


G05.85-0.23 


























1 


17:59:41.27 


-24:03:25.8 










12.13±0.06 


10.88±0.04 


9.08±0.02 


3.58±0.02 




EP 


2 


17:59:49.14 


-24:03:50.6 


14.70±0.04 


11.27±0.03 


9.30±0.02 


7.27±0.01 


6.63±0.01 


6.03±0.01 


5.60±0.01 




2.806 


CII 


3 


17:59:49.88 


-24:03:44.9 




13.97±0.06 


12.91±0.05 


11.87±0.03 


11.48±0.04 


11. 41 ±0.05 






0.963 


CII 


4 


17:59:51.83 


-24:02:04.2 


15.41±0.07 


11.62±0.03 


9.63±0.02 


8.37±0.01 


7.98±0.01 


7.43±0.01 


7.30±0.01 


5.99±0.06 


3.348 


CII 


5 


17:59:47.68 


-24:01:33.0 






13.16±0.06 


11.71±0.02 


11.02±0.01 


10.52±0.01 


10.09±0.05 






CII 


6 


17:59:35.96 


-24:00:43.8 


14.23±0.05 


10.82±0.03 


8.86±0.02 


6.89±0.01 


6.49±0.01 


6.01 ±0.01 


5.35±0.01 




2.777 


CII 


7 


17:59:55.30 


-24:00:39.3 


15.44±0.06 


14.45±0.05 


13.63±0.05 


12.54±0.04 


12.31±0.04 


12.28±0.12 








CII 


8 


17:59:42.43 


-24:00:29.5 


10.28±0.02 


7.61±0.03 


6.28±0.02 


5.68±0.01 


5.50±0.01 


5.14±0.01 


4.96±0.01 


3.40±0.01 


2.176 


cn 


9 


17:59:46.05 


-24:00:15.1 






12.73±0.04 


10.25±0.02 


9.56±0.01 


8.72±0.01 


7.74±0.02 


4.09±0.03 




cn 


10 


17:59:48.08 


-24:00:12.6 


15.30±0.07 


11.35±0.03 


9.23±0.02 


7.28±0.01 


6.63±0.01 


6.15±0.01 


5.57±0.01 


4.07±0.03 


3.482 


cn 


11 


18:00:02.83 


-24:00:07.7 








12.65±0.05 


12.33±0.04 


11.90±0.11 


ll.01iO.18 






CII 


12 


17:59:54.07 


-23:59:42.9 




14.39±0.09 


11.48±0.03 


8.89±0.01 


8.48±0.01 


7.94±0.01 


7.76±0.02 


5.85±0.09 


4.924 


cn 


13 


18:00:02.08 


-23:59:41.2 


15.74±0.09 


14.26±0.06 


13.29±0.06 


12.23±0.03 


11.91±0.04 


11.67±0.10 


11.00±0.28 




0.627 


cn 


14 


18:00:02.14 


-23:59:34.7 






13.83±0.05 


12.38±0.03 


12.07±0.03 


11.49±0.07 


10.59±0.12 






cn 


15 


17:59:58.51 


-23:59:25.3 


15.52±0.06 


14.42±0.05 


13.16±0.04 


12.01±0.03 


11.79±0.04 


11.56±0.07 


11.75±0.32 






cn 


16 


17:59:57.86 


-23:59:12.0 


15.35±0.08 


11.26±0.03 


9.17±0.03 


7.54±0.01 


7.32±0.01 


6.87±0.01 


6.74±0.01 


4.61±0.04 


3.673 


CII 


17 


18:00:02.31 


-23:58:56.1 


15. 85 ±0.09 


12.86±0.04 


11.43±0.03 


10.40±0.01 


10.22±0.01 


9.76±0.01 


9.50±0.03 




2.509 


cn 


18 


17:59:39.24 


-23:58:31.8 






11.99±0.05 


9.77±0.01 


8.95±0.01 


o.zy±u.ui 


/.4D±U.U1 






^TT 


19 


18:00:04.04 


-23:58:04.8 


14.58±0.03 


13.78±0.05 


I3.35±0.06 


12.72±0.05 


12.54±0.04 


12.57±0.14 






0.059 


cn 


20 


17:59:51.58 


-23:57:42.7 


15.20±0.06 


13.88±0.11 


12.92±0.09 


11.73±0.04 


11.33±0.03 


11.18±0.05 


10.63±0.07 




0.400 


cn 


21 


17:59:53.59 


-23:57:40.0 


14.56±0.03 


13.67±0.05 


13.02±0.07 


12.25±0.06 


12.02±0.05 


12.20±0.20 






0.001 ^■ 


cn 


22 


17:59:50.39 


-23:56:59.8 


14.07±0.07 


11.13±0.05 


9.23±0.03 


7.55±0.01 


7.30±0.01 


6.71 ±0.01 


6.17±0.01 




2.102 n 


cn 


23 


17:59:46.44 


-23:56:53.6 


15.26±0.05 


11.49±0.02 


9.45±0.02 


7.67±0.01 


7.27±0.01 


6.71±0.01 


6.30±0.01 




3.278 KH 


cn 


24 


17:59:52.22 


-23:59:03.5 




13.83±0.05 


12.00±0.04 


10.51±0.01 


10.13±0.01 


9.77±0.02 


9.66±0.06 


6.90±0.13 


2.849 S 
Q 


TD 


Note. — Table 2 is published in its entirety in the electronic edition of the Astrophysical Journal. An exceipt is included here to guide the reader in its content and format. 
^ CI=class I protostar, Cn=class II pre-main sequence star, TD=transition disk, EP=embedded protostar 
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TABLES 
YSO Summary 



IRDC 


Class I 


Class n 


Transition 


Embedded 


Total 




Protostars 


PMS stars 


Disks 


Objects 


Number 


G05. 85-0.23 





22(2) 


1 





23(2) 


G06.26-0.51 


3(1) 


26(4) 








29(5) 


G09. 16+0.06 


1 


12(1) 








13(1) 


G09.28-0.15 


2(1) 


15(2) 








17(3) 


G09.86-0.04 


5(3) 


21(1) 


3 


2(1) 


31(5) 


G12.50-0.22 


4(1) 


22(1) 


1 


1(1) 


28(3) 


G23.37-0.29 


8 


36(2) 





2(1) 


46(3) 


G23.48-0.53 


5(4) 


16 








21(4) 


G24.05-0.22 





24(4) 





1 


25(4) 


G34.74-0.12 


4(1) 


28(4) 


2(1) 


1 


35(6) 


G37.44+0.14 


5(1) 


33(1) 


2(1) 





40(3) 


Total 


37(12) 


255(22) 


9(2) 


7(3) 


308(39) 



Note. — The number of objects classified as YSOs for each IRDC field. In parentheses, we indicate the number from each classification that are 
associated with a "clump," which are determined in Section 4.2 and tabulated in Table 5. 



TABLE 4 

SprrzER-mENTiFiED Embedded Protostars: Flux and LuMmosrrY EsTiMArES 



IRDC 


Index 


a 


6 


3.6//m 


4.5jum 


5.8jum 


8jum 


24jum 






number^ 


(J2000) 


(J2000) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(Lo) 


G09.86-0.04 


6 


18:07:36.99 


^20:26:03.9 










7.00±2.47 


>0.05 




7 


18:07:42.12 


-20:23:34.3 










43.64±5.45 


>0.33 


G12.50-0.22 


5 


18:13:41.71 


-18:12:29.6 




0.02± 0.01 






42.93±12.81 


>2.1 


G23.37-0.29 


9 


18:34:54.12 


-08:38:25.5 










33.94± 8.18 


>1.0 




10 


18:35:00.04 


-08:36:57.4 


0.02± 0.01 








18.15±4.56 


>1.5 


G24.05-0.22 


1 


18:35:54.73 


-08:01:30.2 










5.88±1.27 


>0.2 


G34.74-0.12 


5 


18:55:05.20 


-1-01:34:36.2 




0.02± 0.01 






36.13±4.84 


>3.3 



^ In Table 2. 



TABLE 5 
clumpfind parameter summary 



IRDC 


Lower t 
Threshold 


At 


G005.85-0.23 


0.27 


0.20 


G006.26-0.51 


0.27 


0.11 


G009. 16+0.06 


0.19 


0.10 


G009.28-0.15 


0.35 


0.06 


G009.86-0.04 


0.32 


0.11 


G012.50-0.22 


0.31 


0.16 


G023.37-0.29 


0.36 


0.11 


G023.48-0.53 


0.39 


0.09 


G024.05-0.22 


0.22 


0.09 


G034.74-0.12 


0.27 


0.07 


G037.44+0.14 


0.29 


0.17 



Spitzer IRDCs 



TABLE 6 

CLUMPFIND ORIGINAL fg/bo EsmumoN run 



IRDC 








^^iuilip ividaa 




V_.IUIIip olAC 


Notes" 






y > 


I, > 










GO5.85-0.23 


















CI 


23 


-23 


348.5 


0.14 


0.11 






C2 


-22 


-75 


342.5 


0.55 


0.09 






C3 


-33 


-50 


320.7 


0.40 


0.09 






C4 


-60 


-83 


305.8 


0.25 


0.10 






C5 


-16 


-54 


299.0 


0.75 


0.08 






C6 


-6 


-61 


211.8 


0.46 


0.08 






C7 


-42 


-81 


178.4 


0.29 


0.08 


fg 




C8 


-29 


-63 


165.3 


0.93 


0.07 




C9 


78 


18 


108.3 


0.17 


0.08 


fg 




CIO 


-21 


-52 


72.4 


0.58 


0.06 




Cll 


-22 


-61 


64.9 


1.09 


0.05 






C12 


-86 


-93 


47.1 


0.23 


0.07 


5-cn 




C13 


15 


44 


32.9 


0.17 


0.06 


12 -cn 




C14 


-105 


-123 


3.3 


0.14 


0.04 






C15 


-100 


-94 


0.6 


0.13 


0.02 




GO6.26-0.5l 


















CI 


-71 


97 


2226.8 


0.26 


0.18 






C2 


-179 


-25 


963.5 


0.15 


0.15 


17 -cn 




C3 


-99 


93 


820.9 


0.15 


0.15 






C4 


-51 


-45 


683.1 


1.58 


0.11 






C5 


37 


-82 


641.2 


0.21 


0.14 






C6 


-91 


78 


331.7 


0.21 


0.11 






C7 


-154 


-64 


255.4 


0.17 


0.11 






C8 


70 


-112 


240.6 


0.13 


0.11 






C9 


-62 


-64 


233.8 


0.55 


0.10 






CIO 


-59 


-35 


221.8 


0.75 


0.08 






Cll 


-186 


-52 


210.7 


0.17 


0.10 






C12 


-21 


110 


176.8 


0.15 


0.10 






C13 


66 


-117 


175.5 


0.13 


0.11 






C14 


47 


-75 


169.1 


0.20 


0.09 






C15 


91 


-80 


161.1 


0.36 


0.10 


1 - CI 




C16 


-67 


-61 


145.0 


0.49 


0.08 






C17 


-70 


-19 


144.6 


0.45 


0.09 






C18 


-58 


-40 


119.1 


0.88 


0.06 






C19 


70 


-77 


112.5 


0.18 


0.09 






C20 


-71 


56 


109.3 


0.19 


0.09 






C21 


-70 


-30 


107.5 


0.52 


0.07 






C22 


-75 


-39 


106.9 


0.47 


0.08 






C23 


-173 


-44 


102.1 


0.14 


0.09 






C24 


106 


-53 


83.8 


0.14 


0.09 


13 -cn 




C25 


-207 


-38 


77.5 


0.11 


0.09 






C26 


-83 


-86 


77.0 


0.13 


0.09 






C27 


61 


-72 


71.8 


0.19 


0.08 






C28 


-133 


83 


65.4 


0.10 


0.09 






C29 


-105 


-88 


59.6 


0.38 


0.08 






C30 


-61 


-84 


56.3 


0.14 


0.08 






C31 


-167 


-62 


56.0 


0.14 


0.08 


14 -cn 




C32 


-90 


-173 


52.7 


0.48 


0.07 






C33 


-110 


-99 


49.4 


0.20 


0.07 






C34 


-125 


83 


42.5 


0.10 


0.08 






C35 


-135 


-75 


41.8 


0.12 


0.08 






C36 


-110 


70 


37.7 


0.12 


0.07 






C37 


-5 


99 


36.9 


0.13 


0.07 






C38 


-27 


117 


36.6 


0.15 


0.07 






C39 


55 


102 


34.0 


0.14 


0.07 


26 -cn 




C40 


103 


-62 


32.3 


0.15 


0.07 






C41 


-102 


86 


29.9 


0.18 


0.06 






C42 


-178 


-60 


29.7 


0.16 


0.07 


14 - cn 




C43 


65 


-75 


19.8 


0.16 


0.06 






C44 


-25 


121 


18.7 


0.13 


0.06 






C45 


-91 


-49 


17.3 


0.15 


0.06 






C46 


-51 


-84 


15.6 


0.11 


0.06 






C47 


-154 


-12 


14.4 


0.13 


0.06 






C48 


46 


-164 


12.7 


0.10 


0.06 






C49 


-4 


89 


11.0 


0.13 


0.05 






C50 


-129 


80 


10.0 


0.13 


0.05 






C51 


-122 


-86 


9.5 


0.12 


0.05 






C52 


82 


-62 


8.8 


0.13 


0.05 






C53 


65 


117 


8.5 


0.11 


0.05 






C54 


103 


-9 


7.8 


0.13 


0.05 






C55 


-175 


-58 


7.6 


0.15 


0.04 






C56 


-159 


-9 


7.3 


0.10 


0.05 






C57 


-106 


72 


5.5 


0.11 


0.05 





Ragan et al. 



TABLE 6 — Continued 



IRDC 




Aa 


AS 


Clump Mass 


max 


Clump Size 


Notes^ 






(") 


(") 


(Mo) 




(pc) 






C58 


-71 


-177 


4.0 


0.18 


0.04 






C59 


53 


-170 


1.4 


0.10 


0.03 






C60 


-131 


-47 


1.3 


0.10 


0.03 






C61 


77 


ill 


1.2 


0.10 


0.03 






C62 


-77 


-3 


1.0 


0.13 


0.03 






C63 


-102 


-63 


1.0 


0.10 


0.03 




G09. 16+0.06 


















CI 


-63 


-80 


3732.1 


0.52 


0.19 






C2 


-25 


-66 


741.2 


0.43 


0.14 






C3 


-50 


-82 


359.0 


0.76 


0.10 






C4 


-180 


-87 


207.6 


0.17 


0.12 






C5 


-36 


-81 


162.5 


0.52 


0.08 






C6 


1 


-97 


153.6 


0.23 


0.10 


4-cn 




C7 


-213 


-96 


118.6 


0.09 


0.10 






C8 


-106 


-61 


69.7 


0.08 


0.09 






C9 


-24 


-83 


68.3 


0.42 


0.07 






CIO 


-170 


17 


64.1 


0.08 


0.09 






Cll 


-20 


-87 


57.8 


0.32 


0.07 






C12 


-164 


-10 


54.2 


0.11 


0.09 






CIS 


-11 


-111 


51.9 


0.22 


0.07 






C14 


-203 


-1 


38.0 


0.11 


0.08 






C15 


-185 


11 


36.8 


0.11 


0.08 






C16 


-191 


10 


33.2 


0.11 


0.08 






C17 


-123 


-75 


26.0 


0.12 


0.07 






C18 


32 


-118 


23.2 


0.19 


0.07 






C19 


-5 


-130 


16.3 


0.16 


0.06 






C20 


-194 


-31 


10.8 


0.08 


0.06 






C21 


-2 


-82 


10.7 


0.19 


0.05 


4-cn 




C22 


-119 


-38 


4.4 


0.08 


0.05 






C23 


-116 


-45 


4.0 


0.09 


0.04 






C24 


87 


-92 


2.9 


0.11 


0.04 






C25 


-100 


-40 


1.7 


0.08 


0.04 






C26 


-69 


-23 


1.5 


0.11 


0.04 






C27 


-104 


-39 


1.3 


0.08 


0.03 




GO9.i8-0.i5 


















CI 


-77 


1 


1036.9 


0.55 


0.13 


fg 




C2 


-64 


-21 


636.1 


0.77 


0.10 


fg 




C3 


-59 


-14 


492.4 


0.61 


0.11 




C4 


-55 


-43 


343.3 


0.71 


0.10 






C5 


-81 


-12 


339.3 


0.50 


0.10 






C6 


-49 


-56 


283.5 


0.49 


0.09 


fg 




C7 


-50 


-32 


238.3 


0.71 


0.09 




C8 


-89 


36 


234.6 


0.21 


0.11 


14 -cn 




C9 


-127 


-25 


169.2 


0.31 


0.10 






CIO 


-80 


23 


119.8 


0.28 


0.09 






Cll 


-36 


-125 


119.3 


0.29 


0.09 


fg 




C12 


-55 


-98 


118.8 


0.23 


0.09 




C13 


-54 


-38 


115.2 


0.67 


0.07 






C14 


-37 


-112 


110.2 


0.27 


0.09 


fg 




C15 


-51 


-70 


106.8 


0.40 


0.08 




C16 


-43 


-39 


96.3 


0.41 


0.08 






C17 


-59 


-76 


87.0 


0.24 


0.08 






C18 


-34 


-98 


73.0 


0.30 


0.08 






C19 


-16 


-147 


72.5 


0.20 


0.08 






C20 


-139 


-25 


70.0 


0.23 


0.08 


fg 




C21 


-38 


-104 


65.5 


0.29 


0.07 


fg 




C22 


8 


-34 


60.8 


0.20 


0.08 




C23 


-101 


-9 


60.2 


0.22 


0.08 






C24 


-97 


-20 


59.0 


0.21 


0.07 






C25 


-102 


-23 


58.2 


0.19 


0.08 


fg 




C26 


-96 


27 


51.5 


0.27 


0.07 




C27 


-34 


-137 


48.4 


0.20 


0.08 


10 -cn 




C28 


-37 


-83 


46.8 


0.38 


0.07 






C29 


-25 


-65 


45.8 


0.23 


0.07 






C30 


-57 


-92 


45.3 


0.21 


0.07 






C31 


-40 


-90 


44.2 


0.38 


0.06 






C32 


-80 


27 


43.6 


0.29 


0.07 






C33 


4 


-26 


41.5 


0.19 


0.07 


2 -CI 




C34 


-76 


-33 


37.6 


0.16 


0.07 






C35 


-26 


-45 


31.7 


0.20 


0.07 


fg 




C36 


-41 


-69 


31.6 


0.28 


0.06 




C37 


-32 


-74 


30.9 


0.23 


0.06 






C38 


-36 


-61 


28.1 


0.28 


0.06 






C39 


-102 


-15 


27.1 


0.22 


0.06 






C40 


-53 


-80 


26.2 


0.26 


0.06 






C41 


-49 


-49 


25.8 


0.41 


0.05 





Spitzer IRDCs 
TABLE 6 — Continued 



IRDC 




Aa 
(") 


AS 
(") 


Clump Mass 

(Mo) 


'^max 


Clump Size 

(pc) 


Notes" 




C42 


-32 


-J J 


1 Q 1 


ft 99 
u.zz 


06 






C43 




-oo 


1 R 1 

io. i 


ft 9n 

D.ZD 


\j.\jyj 






C44 


1 Q 




i J.D 


ft 9S 


05 


fg 




C45 


-20 


-OD 


i J.D 


23 


05 






AO 


1 

-ID 


1 7 

ID. / 


ft 1 1 


OS 






PAT 


- IKJ 




ID. J 


ft 90 

D.ZU 


U.UJ 






C48 


-107 


-50 


10.0 


0.14 


0.05 


fu 




C49 


-17 


-161 


9.6 


0.16 


0.05 




C50 


-4 


-146 


6.7 


0.16 


0.05 






C51 


-12 


31 


4.6 


0.13 


0.04 






C52 


-58 


7 


4.5 


0.11 


0.04 






C53 


-110 


-39 


4.2 


0.15 


0.04 






C54 


-12 


-93 


2.6 


0.14 


0.04 




C55 


-139 


12 


2.6 


16 


04 






C56 


-141 


44 


2.0 


0.12 


0.04 






C57 


-32 


-157 


1.8 


0.13 


0.03 




C58 


-113 


45 


1.8 


0.14 


0.03 
























1 < 

- LJ 


-Do 


9QQ ft 


42 


08 






r^9 


-DZ 


S9 
-oZ 


1 S'^ 7 
1 Oj. / 


ft 1 f\ 


n OR 

\).\)o 






y^j 


-37 


-71 


174 3 


20 


07 


Ig 




C4 


-101 


S7 
-J / 


ID / .7 


50 


07 


Ig 




C5 


-74 


-44 


1 Q 

i JD.7 


1 57 


06 






Q9 
-7Z 


-42 


19ft '^ 

iZD. J 


ft AS 


C\f\ 
D.DD 






Cl 


^9 

-JZ 


-Jo 


1 17 A 


ft 97 
D.Z / 


Cif, 

u.uo 






C9. 
^_-o 


-42 


-Jo 


1 Ifi 

1 ID. J 


41 


O^i 
u.uo 


6 - EP 






Q 


fn 

-D / 


119^ 

i IZ. J 


ft ^9 

D. jZ 


07 

U.U / 


Ig 






^'^ 


-44 


IHQ Q 


42 


Ofi 
u.uo 




PI 1 




-jj 


09 4 


ft '^9 

D.JZ 


0*^ 

U.UJ 






C12 


-85 


-61 


82.7 


0.30 


0.06 






C13 


-126 


-28 


80.4 


0.59 


0.06 






Pii 


-OJ 


-39 


Rft ft 

OD.D 


ft 77 


OS 

U.UJ 






PI S 


-24 


-1 1 


/ D.-J 


18 


0^^ 






PI 


-1 14 


-106 


70 6 


18 


06 






PI 7 


-72 


-27 


65 9 


26 


06 






C18 


-1 15 


-38 


63 9 


64 


05 






C19 


-15 


-40 


56.1 


17 


06 






C20 


9/^ 
-zo 


-Di 


ilR S 
to. J 


ft "^R 

U.-JO 


05 






P9 1 


-J 


S7 
-o / 


il^ 9 


ft 1 1 


0^1 






P99 


1 < 

- i J 


-JD 


40 8 


22 


OS 


Ig 




P9'^ 


-1 14 


19^ 

- iZj 


OD.f 


ft 17 
D. i / 


OS 

U.UJ 




C24 


-i-jy 


7 




24 


OS 

U.UJ 


10 PTT 




P9'^ 


-1 J / 


-121 


97 7 


ft 9^^ 

D.ZJ 


OS 

U.UJ 


fo 

Ig 




P9fi 




-122 


9fi 

ZD.D 


22 


OS 

U.UJ 




C27 


-121 


-31 


23.8 


0.56 


0.04 






C28 


27 


-65 


21.9 


0.32 


0.05 


fo 

Ig 




C29 


_5 


-39 


19.1 


0.15 


0.04 


5 - CI 




C30 


-52 


-45 


18.0 


0.28 


0.04 






C31 


-45 


-52 


17 8 


26 


0.04 






C32 


-184 


-125 


12.9 


17 


0.04 






C33 


-160 


-131 


6.5 


0.13 


04 






C34 


-157 


_1 


4.7 


0.12 


0.03 






C35 


-3 


-116 


7 


12 


02 






C36 


-201 


j-f 


4 


1 1 


02 






















PI 


-70 


-30 


1418 6 


39 


14 






P9 

L-Z 


-J i 


-61 


1293 3 


44 


14 


2 - CI 




C3 


-41 


-59 


838 6 


62 


0.12 










-41 


488 8 


98 


09 


J - Dr, LA - V^ll 




P'^ 


1 Qf\ 




JO J.t 


ft Af\ 
D.tD 


Ift 
U. ID 






C6 


-104 


-31 


373.2 


0.23 


0.11 


fe 




C7 


-51 


-50 


333.5 


0.95 


0.08 


2 - CI 




C8 


-78 


-12 


204.6 


0.40 


0.09 


14 - cn 




C9 


-42 


-43 


179.0 


1.16 


0.07 


5 - FP 1 2 - rn 




CIO 


-190 


49 


165.1 


0.32 


0.08 






Cll 


-131 


-47 


99.9 


0.26 


0.08 






C12 


80 


-146 


36.2 


0.21 


0.07 


fg 




C13 


-32 


-43 


31.1 


0.25 


0.06 




C14 


-32 


-112 


21.9 


0.22 


0.06 






C15 


-35 


-101 


3.6 


0.20 


0.04 




G23.37-0.29 


















CI 


45 


43 


5199.6 


1.64 


0.18 






C2 


45 


37 


3143.6 


1.28 


0.17 


9-EP 




C3 


-18 


183 


417.3 


0.32 


0.13 






C4 


-44 


139 


317.3 


0.82 


0.10 





Ragan et al. 
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IRDC 




Aor 


A5 


Clump Mass 




Clump Size 


Notes^ 






D 


(") 


(Mo) 




(pc) 








- it) 


1Q7 


231 2 


58 


Oil 








-44 




1 RS 

100.7 


28 


10 
0. 10 






n 


o / 


24 


1 9 
100. z 


9 1 
o.z i 


oil 






C8 


48 


70 


167 5 


24 


10 








7Q 


1 (17 


1 9R 

iZo.D 


1 8 
0. 10 


10 
0. 10 








l7Z 


900 


197 S 

IZ / ,J 


1 10 

1.17 


00 






Cll 


-50 


167 


94.0 


0.26 


0.09 






C12 


-10 


101 


81 4 


'^1 

U. Jl 


08 
0.00 




C13 


39 


112 


78.8 


0.27 


0.08 








7 
1 


-00 


78 9 
/ o.z 


x\ 


08 

U.Uo 






C15 


-26 


206 


67.3 


0.17 


0.09 






V- ID 


u 


1 SI 


OJ.M- 


20 


09 


10 PTT 




pn 




194 


^^1 
Oi .u 


18 


08 






PI 8 

V- I o 


124 


1 '^ 


10 8 


14 


07 






PI Q 


- 


126 


90 1 
zy. J 


21 


07 






p?n 


171 




28 


98 
O.ZO 


c\f 
0.00 






P71 


1 n 


1 OA 


97 7 


9A 


07 

u.u / 






P79 


-52 


186 


24 6 


98 

O.Zo 


06 


fa 








1 09 

lUZ 


9A 9 
Z'+.Z 


91 
O.Z 1 


U.UO 




P9A 


«: 

J 




99 9 
ZZ.Z 


91 
O.Z 1 


0^^ 

U.UO 






P9S 




1 99 

IZZ 


9 1 1 
Zi . i 


99 
U.ZZ 


U.UO 




v^zu 


- /o 


1 1 8 


90 1 
ZO. J 


17 
0. 1 / 


c\f 

U.UO 






v^z / 


1 Q9 

i7Z 


202 


1 8 

10. D 


1 40 


OS 
u.u J 






P9R 


n 
u 


-Of 


18 2 


90 


dfx 
U.UO 




P9Q 


-J7 


1 1 S 


1 S 


21 


Ofi 

U.UO 






C30 


-47 


199 


15.1 


0.18 


0.06 






C31 


74 


-65 


14.9 


0.16 


0.06 






P'^9 

V^JZ 


^4 
-jt 


211 


1 A 7 

It. / 


1 8 

0. 10 


Ofi 

U.UD 






C33 


49 


98 


13.9 


0.25 


0.05 


9 - EP 




C34 


53 


88 


13 7 


20 


06 






C35 


15 


114 


12 4 


18 


06 


fg 






inn 


172 


12 1 


1 f, 
0. 10 


c\f 

U.UO 




CXI 


zu 


114 


12 


1 s 

0. 1 J 


Of 

U.UO 


fa 




P'^8 


1 Qfs 


i jj 


10 fx 
lU.O 


1 s 

0. 1 J 


Of 

U.UO 






JO 


70 


8 9 
o.z 


17 


05 






C40 


71 


-43 


8.1 


19 


05 






C41 


-62 


188 


6.2 


0.17 


0.05 






C42 


-63 


199 


6.0 


0.16 


0.05 






C43 


164 


140 


5.4 


0.16 


0.05 






C44 


58 


86 


5.4 


13 


05 


Ig 




C45 


-83 


185 


4.3 


15 


04 




C46 


-109 


70 


4.2 


0.14 


0.04 






C47 


sn 


-58 


4.0 


0.15 


0.04 






C48 


-29 


38 


3.8 


0.18 


0.04 


33 - cn 




C49 


-11 


81 


2.6 


0.13 


0.04 


Ig 




C50 


-56 


202 


1.8 


0.14 


0.04 


KJZ/J.HO — U.J J 


















PI 


78 


26 


1274 2 


82 


13 








68 


32 


775 


71 


12 






P'^ 


104 


29 


436 3 


33 


12 






CA 


JZ 


37 


JOJ.O 


76 


00 
u.uy 






P'> 


tj 




249 7 


91 


08 


A PT 

f - V^l 




p/^ 


SO 


s^^ 


247 5 


58 


09 






C7 


_43 


AS 


206 2 


41 


09 






C8 


33 


40 


204 1 


72 


08 


A PT 




PQ 


12 


42 


202 3 


58 


09 






Pin 


60 


5 1 


186 1 


37 


09 






CI 1 


-79 


-72 


182 4 


31 


09 


1 - PT 9 - PT 1 - PT 




C12 


1Q 


SO 

J7 


178 
1 / o.u 


AS 


00 

U.U7 






PT^ 


-DZ 


^i9 

-DZ 


154 1 




00 

U.U7 






Pll 


Q1 


n 

D 


142 7 


90 


00 

U.U7 






PI 




S7 
-J / 


190 fx 


AO 


08 

U.Uo 












190 


11 

O.J 1 


00 

U.U7 






P17 


9'^ 

ZJ 


<1 

Jl 


101 2 


fil 

O.Dl 


(\f 

U.UD 






PIS 


-12 


-24 


Ofi 7 


98 

O.ZO 


08 

U.Uo 






PIQ 




so 


0^ f\ 


Af, 

O.f D 


07 

U.U / 






C20 


Q 
o 


7S 


y 


28 


08 






C21 


44 


44 


83.3 


0.59 


0.06 


4-CI 




C22 


-24 


-29 


79.2 


0.26 


0.08 






C23 


-54 


-64 


69.7 


0.38 


0.07 






C24 


116 


30 


64.0 


0.19 


0.08 






C25 


-1 


57 


47.5 


0.26 


0.07 






C26 


28 


54 


46.4 


0.52 


0.05 






C27 


4 


37 


40.4 


0.46 


0.06 






C28 


-91 


101 


39.0 


0.23 


0.07 






C29 


-9 


-12 


38.6 


0.19 


0.07 





Spitzer IRDCs 
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IRDC 




Aa 


A5 
f") 


Clump Mass 

V-'-'-'O/ 


'max 


Clump Size 


Notes" 




C30 


105 


-10 


38.0 


0.29 


0.07 






C31 


36 


40 


27.1 


0.72 


0.04 


4-CI 




C32 


68 


68 


18.0 


0.25 


0.06 






C33 


123 


42 


14.5 


0.17 


0.05 






C34 


52 


80 


12.9 


0.17 


0.05 






C35 


59 


70 


9.6 


0.21 


0.05 


fg 




C36 


-59 


-44 


5.0 


0.16 


0.04 


fg 




C37 


66 


3 


4.5 


0.17 


0.04 




C38 


-21 


89 


2.4 


0.17 


0.03 




G24.05-0.22 


















CI 


44 


55 


598.9 


0.70 


0.12 






C2 


32 


37 


411.5 


0.58 


0.11 


7-cn 




C3 


26 


56 


363.0 


0.37 


0.12 






C4 


31 


83 


343.8 


0.29 


0.13 


fg 




C5 


44 


20 


277.8 


0.28 


0.12 




C6 


40 


84 


217.8 


0.21 


0.12 


fg 




C7 


36 


183 


192.1 


0.10 


0.13 


13 - CII 




C8 


38 


46 


184.5 


1.10 


0.08 


10 -cn 




C9 


47 


128 


181.7 


0.14 


0.12 






CIO 


40 


165 


179.6 


0.14 


0.12 






Cll 


173 


181 


131.7 


0.17 


0.11 






C12 


52 


42 


122.2 


0.36 


0.09 






C13 


56 


119 


118.1 


0.19 


0.10 






C14 


45 


73 


99.4 


0.23 


0.09 


fg 




C15 


69 


124 


71.9 


0.16 


0.09 




C16 


28 


160 


70.7 


0.11 


0.10 






C17 


66 


113 


61.5 


0.20 


0.09 






C18 


44 


45 


54.0 


0.44 


0.06 






C19 


144 


205 


53.6 


0.11 


0.09 






C20 


59 


37 


47.4 


0.23 


0.08 






C21 


130 


208 


45.2 


0.13 


0.09 


fg 




C22 


152 


173 


33.5 


0.12 


0.08 




C23 


136 


177 


27.6 


0.11 


0.08 






C24 


28 


45 


27.1 


0.31 


0.06 






C25 


-90 


173 


23.9 


0.19 


0.07 


fg 




C26 


128 


186 


19.3 


0.11 


0.07 


fg 




C27 


184 


172 


15.6 


0.11 


0.07 




C28 


24 


127 


15.1 


0.13 


0.07 


8 -Cll 




C29 


71 


59 


14.6 


0.13 


0.07 


8-cn 




C30 


191 


179 


14.1 


0.10 


0.07 






C31 


68 


32 


12.8 


0.19 


0.06 






C32 


148 


175 


8.9 


0.10 


0.06 






C33 


-83 


-24 


7.8 


0.14 


0.05 






C34 


-94 


173 


4.8 


0.19 


0.05 


fg 




C35 


-100 


167 


4.4 


0.14 


0.05 


fg 




C36 


76 


37 


3.8 


0.14 


0.05 




C37 


-106 


176 


3.7 


0.11 


0.05 


fg 




C38 


-75 


-42 


2.0 


0.10 


0.04 




C39 


7 


184 


1.8 


0.09 


0.04 




G34.74-0.12 


















CI 


-62 


-26 


1411.1 


0.66 


0.14 






C2 


-59 


-49 


353.2 


0.42 


0.11 






C3 


-58 


-39 


309.7 


0.42 


0.10 






C4 


-5 


-83 


282.1 


0.19 


0.13 


14 - CII 




C5 


-23 


-60 


223.0 


0.49 


0.10 






C6 


-11 


-62 


185.3 


0.32 


0.10 


14 - Cll 




C7 





-108 


157.0 


0.09 


0.12 






C8 


-36 


-54 


150.8 


0.35 


0.10 


35 - TD 




C9 


-39 


-40 


147.7 


0.50 


0.08 






CIO 


-19 


-40 


146.8 


0.35 


0.09 






Cll 


-26 


-76 


144.9 


0.26 


0.10 






C12 


-59 


-10 


134.7 


0.22 


0.10 






C13 


-37 


-17 


113.7 


0.22 


0.10 






C14 


-47 


-26 


92.3 


0.29 


0.08 






C15 


-89 


-24 


89.8 


0.14 


0.10 


fg 




C16 


-39 


-83 


86.6 


0.12 


0.10 




C17 


-44 


-5 


83.4 


0.21 


0.09 


fg 




C18 


-32 


-39 


74.7 


0.29 


0.08 




C19 


-42 


-28 


69.4 


0.27 


0.08 






C20 


7 


-59 


63.1 


0.16 


0.09 






C21 


-6 


-38 


57.6 


0.42 


0.07 


3 -CI 




C22 


-15 


-89 


54.6 


0.15 


0.09 


fg 




C23 


-122 





54.4 


0.19 


0.09 


18 -cn 




C24 


-26 


-48 


54.3 


0.34 


0.07 






C25 


-108 


-15 


52.5 


0.17 


0.09 






C26 


-48 


-15 


50.2 


0.22 


0.07 
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IRDC 


Aa 


AS 


Clump Mass 


'max 


Clump Size 


Notes^ 




V / 


(") 










C27 


-28 


-106 


46.2 


0.11 


0.09 


fg 


C28 


-8 


-43 


42.6 


0.38 


0.06 


3 - CI, 16 - Cll 


C29 


-5 


-55 


41.7 


0.26 


0.07 




C30 


-36 


-44 


41.0 


0.50 


0.06 




C31 


-40 


5 


40.5 


0.12 


0.08 




C32 


-8 


-33 


34.3 


0.22 


0.07 




C33 


-28 


-91 


34.0 


0.15 


0.08 




C34 


-14 


-49 


32.8 


0.30 


0.06 


16 -cn 


C35 


18 


-73 


30.5 


0.09 


0.08 




C36 


-81 


-114 


29.8 


0.13 


0.08 




C37 


-141 


19 


28.6 


0.17 


0.07 


22 -cn 


C38 


-29 


-113 


28.3 


0.09 


0.08 




C39 


1 


-53 


28.0 


0.20 


0.07 




C40 


12 


-103 


25.7 


0.14 


0.07 




C41 


-12 


-83 


25.2 


0.14 


0.07 




C42 


-84 


-106 


23.5 


0.16 


0.07 


fg 


C43 


-42 


63 


23.4 


0.15 


0.07 


fg 


C44 


-127 


-60 


18.1 


0.16 


0.07 


C45 


-56 


10 


16.7 


0.13 


0.07 




C46 


-92 


-16 


15.4 


0.12 


0.06 


fg 


C47 


-60 


-62 


13.8 


0.14 


0.06 


C48 


-47 


-94 


13.0 


0.11 


0.06 




C49 


-130 


9 


12.4 


0.11 


0.06 




C50 


-21 


-86 


11.4 


0.13 


0.06 




C51 


6 


-113 


8.7 


0.10 


0.06 




C52 


-48 


57 


8.6 


0.13 


0.06 




C53 


-25 


-92 


5.8 


0.14 


0.05 


fg 


C54 


-131 


-81 


2.8 


0.09 


0.04 


C55 


-110 


-5 


2.1 


0.09 


0.04 




G37.44+0.14 














CI 


52 


159 


352.7 


0.19 


0.10 




C2 


64 


48 


258.3 


0.45 


0.08 


2-cn 


C3 


77 


43 


218.6 


0.40 


0.08 




C4 


-86 


36 


203.7 


0.30 


0.08 




C5 


45 


141 


170.1 


0.17 


0.08 




C6 


70 


196 


130.4 


0.17 


0.08 


fg 


C7 


35 


206 


119.3 


0.14 


0.08 


C8 


61 


211 


117.1 


0.29 


0.07 


fg 


C9 


30 


47 


116.4 


0.80 


0.06 


CIO 


-101 


23 


110.0 


0.35 


0.07 




Cll 


18 


203 


95.7 


0.16 


0.07 




C12 


24 


58 


89.8 


0.92 


0.06 


40 -TD 


C13 


51 


56 


88.6 


0.59 


0.06 




C14 


47 


205 


68.6 


0.23 


0.06 




C15 


52 


210 


60.4 


0.25 


0.06 




C16 


-111 


31 


53.4 


0.81 


0.05 




C17 


-102 


38 


52.7 


0.38 


0.05 




C18 


51 


72 


44.8 


0.59 


0.05 




C19 


44 


16 


35.4 


0.15 


0.06 




C20 


5 


16 


34.6 


0.32 


0.05 




C21 


42 


56 


26.8 


0.60 


0.04 




C22 


-114 


40 


25.5 


0.98 


0.04 




C23 


17 


42 


24.7 


0.23 


0.05 




C24 


41 


77 


22.5 


0.64 


0.04 




C25 


91 


39 


20.9 


0.18 


0.05 




C26 


33 


131 


17.8 


0.13 


0.05 




C27 


28 


14 


17.5 


0.16 


0.05 




C28 


97 


33 


16.2 


0.42 


0.04 




C29 


42 


23 


15.8 


0.13 


0.05 


20 -cn 


C30 


13 


48 


15.3 


0.13 


0.05 




C31 


36 


61 


14.2 


0.60 


0.03 




C32 


105 


34 


13.3 


0.32 


0.04 




C33 


31 


60 


13.0 


0.87 


0.03 




C34 


34 


3 


12.6 


0.14 


0.04 




C35 


11 


18 


11.3 


0.23 


0.04 




C36 


37 


55 


10.3 


0.55 


0.03 




C37 


5 


34 


10.0 


0.13 


0.04 




C38 


112 


33 


10.0 


0.28 


0.04 




C39 


53 


82 


8.7 


0.21 


0.04 




C40 


33 


105 


7.2 


0.13 


0.04 


fg 


C41 


29 


29 


6.7 


0.14 


0.04 


C42 


-52 


88 


5.9 


0.83 


0.03 




C43 


37 


20 


5.0 


0.15 


0.03 




C44 


-108 


-19 


4.4 


0.12 


0.03 


fg 


C45 


68 


66 


3.7 


0.12 


0.03 


C46 


-112 


37 


2.4 


0.96 


0.02 
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IRDC 


Aa 


A5 


Clump Mass 


^max 


Clump Size 


Notes^ 




(") 


V ) 










C47 


85 


54 


1.9 


0.12 


0.03 




C48 


8 


37 


1.8 


0.13 


0.03 




C49 


-111 


134 


1.7 


0.26 


0.03 




C50 


-46 


82 


1.7 


0.27 


0.02 




C51 


-59 


88 


1.1 


0.19 


0.02 




C52 


110 


39 


1.0 


0.30 


0.02 




C53 


136 


-6 


0.9 


0.15 


0.02 




C54 


-103 


121 


0.9 


0.21 


0.02 


fg 


C55 


90 


16 


0.5 


0.11 


0.02 



^ Numbers indicate which stars from Table 2 are associated with a given clump. The YSO type (CI = Class I; Cn = Class H; EP = Embedded Protostar; TD = Transition Disk) is also 
listed. The "fg" denotation indicates that a foreground (or background) star in the field may contaminate the properties listed for that clump. 



